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Figure 
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26. 

Differential  cross  sections  for  Born 
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(described  in  text,  pp.  IO7-IO8). 
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THE  ELASTIC  SCATTERING  OP  POSITIVE  AND 
NEGATIVE  PIONS  BY  NUCLEI 

I ; Introduction; 

An  investigation  of  the  elastic  scattering  of  positive 
and  negative  pions  by  various  nuclei  is  now  in  progress.  The 
elastic  scattering  of  60  Mev  positive  and  negative  pions  has 
been  previously  investigated  by  Byfield  et  al.,^  using  carbon 
plates  in  an  expansion  cloud  chamber  placed  in  the  Nevis  60  Mev 
meson  beams.  An  analysis  of  the  difference  of  the  positive 
and  negative  pion  elastic  scattering  using  a modified  optical 
model  and  taking  into  account  the  interference  of  the  Coulomb 
potential  with  the  nuclear  potential  in  a Born  approximation 
led  to  the  conclusion  that  the  nuclear  potential  experienced 
by  the  Incident  meson  was  most  likely  attractive. 

p 

An  independent  analysis  by  Peaslee  ' which  essentially 
represented  an  addition  of  the  individual  nucleon  scattering 
contributions  came  to  the  same  conclusion  and  therefore  im- 
plied that  the  sign  of  the  p wave  phase  shift  (which  seems  to 
be  mostly  responsible  for  the  low  energy  meson  scattering 
phenomena)  was  positive. 

A later  cloud  chamber  investigation  of  the  elastic 
scattering  of  12 $ Mev  negative  pions  by  carbon  and  lead  per- 

"i 

formed  by  Kessler  and  Lederman-'  was  also  consistent  with  an 
optical  model  analysis  using  an  attractive  nuclear  potential. 
However,  the  lack  of  corresponding  positive  pion  data  did  not 


2 


allow  as  sensitive  a determination  of  the  sign  of  the  poten- 
tial . 

Similar  cloud  chamber  experiments  for  negative  mesons 

on  carbon  were  carried  out  by  Shapiro^-  and  analyzed  by 

t 6 

Bethe  and  Wilson.  Tracy,  with  poor  statistics  has  measured 
the  elastic  scattering  of  x+  - mesons  on  aluminum  by  using  a 
cloud  chamber  while  Shutt  et  al,"^  using  a diffusion  chamber, 
have  measured  the  scattering  of  pions  on  helium. 

O 

Isaacs  et  al°  also  investigated  the  elastic  scattering 
of  60  Mev  negative  pions  from  carbon  using  scintillation 
counter  techniques  and  obtained  much  better  statistics.  The 
data  was  also  consistent  with  a negative  well  but  here  again 
due  to  the  lack  of  positive  data  the  determination  of  the  sign 
of  the  potential  was  not  very  sensitive. 

Following  this  early  work  a systematic  survey  of  the 
elastic  scattering  of  both  positive  and  negative  pions  of 
various  energies  for  various  elements  (both  heavier  and 
lighter  than  carbon)  was  undertaken  at  Nevis.  Scintillation 
counter  techniques  were  employed  and  make  possible  angular 
resolutions  and  statistical  accuracies  far  superior  to  those 
obtained  in  the  cloud  chamber  experiments.  This  paper  re- 
ports on  the  initial  result  in  this  program  for  80  Mev  n*  and  x 
mesons  on  aluminum. 

The  purpose  of  these  experiments  was: 

(l)  To  investigate,  with  improved  statistical  accuracy 
and  angular  resolution,  the  angular  distribution  of  x mesons 
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from  nuclei  as  an  important  phenomenon  in  high  energy  physics. 

In  particular,  it  is  of  interest  to  investigate  as  carefully 
as  possible  the  differences  between  x+  and  it”  scattering  in 
nuclei  with  roughly  equal  numbers  of  neutrons  and  protons  to 
observe  interference  effects  between  "meson  forces'1  and 
"electromagnetic  forces"  in  the  scattering. 

(2)  There  is  some  interest  in  seeing  how  well  the  scatter- 
ing can  be  matched  by  an  "optical  model",  where  the  effect  of 
nuclear  interactions  is  represented  by  a complex  potential 
inside  the  nucleus.  Of  most  interest  is  the  question  of  the 
possibility  of  using  a constant  complex  potential  inside  the 
nucleus  and  a coulomb  field  outside.  The  imaginary  part  of 

the  potential  attenuates  the  wave  to  take  into  account  ab- 
sorption, inelastic,  and  other  possible  incoherent  processes, 
and  can  be  alternately  expressed  in  terms  of  a mean  free  path 
in  nuclear  matter  for  such  processes. 

(3)  To  compare  meson  angular  scattering  with  that  of 
other  probe  particles,  such  as  electrons,  neutrons,  and  protons, 
to  give  information  on  the  distribution  of  nuclear  matter  as  a 
function  of  distance  from  the  center  of  the  nucleus. 

It  is  to  be  expected  that  such  scattering  will  be  of  con- 
tinuing interest  to  physics  and  that  the  accuracy  of  such 
measurements  will  be  steadily  improved  to  allow  comparison  with 
more  refined  theoretical  analysis. 

(4)  It-  is  of  interest  to  see  if  the  scattering  can  be 
explained  in  terms  of  the  known  results  for  the  scattering  by 
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single  nucleons  with  the  application  of  the  charge  inde- 
pendence hypothesis.  This  need  not  apply,  since  the  properties 
of  nuclear  matter  may  not  be  able  to  be  represented  in  terms 
of  individual  processes  between  pairs  of  particles  alone. 

PLAN  OF  THE  EXPERIMENT 


II. 

A.  General 

The  Nevis  3 85  Mev  proton  cyclotron  produces  a beam  of 
circulating  protons.  When  a beryllium  target  is  inserted 
inside  the  cyclotron,  protons  striking  the  target  induce 
nuclear  reactions  which  emit,  among  other  products  of  reaction, 
negative  and  positive  mesons  of  varying  momenta.  These  charged 
particles  are  analyzed  and  focussed  by  the  fringing  field  of 
the  cyclotron  magnet  into  beams  of  equal  momentum  particles. 
Thus  with  the  cyclotron  magnetic  field  In  a given  direction, 
particles  of  one  sign  (say  negative  mesons)  emerge  in  the 
direction  of  the  experimental  area.  Channels  are  cut  in  the 
cyclotron  shielding  wall  (8>  of  lead  and  iron)  to  allow  the 
mesons  to  enter  the  experimental  area  (see  Figure  1.)  To 
obtain  the  opposite  sign  meson  particle  the  magnetic  field 
of  the  cyclotron  is  reversed.  Then  positive  mesons  of  the 
same  momenta  as  the  negative  ones  will  be  focussed  through 
the  same  channel  in  the  identical  manner.  Mesons  of  180  Mev/c 
mean  momentum  emerged  from  the  channel  chosen  in  this  experi- 
ment. This  channel  was  chosen  as  it  corresponded  to  the 
highest  energy  channel  from  which  both  positive  and  negative 
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mesons  were  available  with  good  intensities. 

Upon  emerging  from  the  shielding  wall  the  meson  beam 
was  deflected  and  focussed  by  an  instrument  magnet.  The 
beam  was  then  defined  by  a scintillation  counter  telescope 
consisting  of  two  scintillation  detectors,  about  three  feet 
apart,  connected  to  a coincidence  circuit  that  effectively 
recorded  the  number  of  particles  passing  throxigh  both  of  the 
detectors.  Following  this  telescope  and  in  line  with  it,  the 
target  to  be  studied,  1/2"  or  l/l±n  thick  aluminum  5”  by  3" 
in  size,  was  mounted  on  a turret  stand.  Most  of  the  incident 
beam  striking  the  target  passed  through,  but  a small  percentage 
interacted  in  the  aluminum  target.  Part  of  the  interaction 
was  elastic  and  yielded  elastically  scattered  particles  into 
all  angles  from  0 to  180  degrees  with  the  incident  beam 
direction.  To  detect  these  elastically  scattered  particles 
two  more  scintillation  detectors  were  mounted  on  the  arms  of 
a rotatable  scattering  turret  (Figure  2)  that  allowed  the 
detectors  to  be  rotated  in  a vertical  plane,  to  any  angle 
with  the  incident  beam.  By  placing  the  coincidence  output 
of  these  last  two  detectors  in  coincidence  with  the  output 
of  the  monitoring  telescope  it  was  possible  to  measure  the 
number  of  beam  particles  scattered  at  a given  angle.  Sufficient 
absorber  was  placed  between  the  last  two  detectors  so  as  to 
allow  the  elastically  scattered  particles  to  be  predominantly 
detected  and  to  absorb  any  secondary  protons  or  inelastically 


scattered  mesons 
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B.  Focussing  Magnet 

The  focussing  magnet  served  several  purposes.  As  a 
bending  magnet  it  supplemented  the  momentum  selection  of  the 
cyclotron  fringing  field.  It  also  deflected  the  beam  out  of 
line  with  the  channel  opening.  The  opening  views  the  interior 
of  the  cyclotron  enclosure,  and  a rather  high  background  is 
present  along  the  line  of  the  main  channel  opening.  This  back- 
ground consists  of  garnma  rays,  neutrons,  and  some  charged 
secondary  particles  which  are  mainly  not  deflected  the  same 
as  the  meson  beam.  Triangular  shaped  pole  pieces  were  chosen 
to  obtain  an  increase  in  intensity  through  focussing.  The 
magnet  acted  as  a converging  lens  in  the  vertical  direction, 
and  as  a defocus sing  lens  in  the  horizontal  direction.  A 
net  gain  of  intensity  of  a factor  of  three  was  obtained 
through  the  monitoring  telescope.  Raising  or  lowering  a 
l/2"  crystal  at  the  focal  point  by  7/8'*  reduced  the  intensity 
by  nearly  a factor  of  three.  A further  gain  was  then  realized 
by  using  rectangularly  shaped  counters,  similar  to  the  shape 
of  the  source,  with  the  smaller  dimension  vertical.  The 
scattering  was  measured  through  angles  in  a vertical  plane. 

C.  Scintillation  Detectors 

Stilbene  crystals  were  used  as  the  scintillators  with 
a decay  time  of  less  than  8 x 10-<^  seconds.  The  first  crystal 
of  the  telescope,  XI,  was  4"  x 2n  x l/8n • the  second,  X2, 
was  4"  x 1/2”  x 1/8"  in  the  first  runs  and  was  later  changed 
to  3"  x 3/4"  x 3/16" . The  crystals  mounted  on  the  rotating 
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arm  were  both  I4."  x 2";  the  nearer,  X3,  being  l/8”  thick,  and 
the  last,  Xlj.,  being  I/I4."  thick. 

For  maximum  efficiency  in  counting  particles  scattered 
at  a given  angle,  an  annular  counter,  subtending  an  angle 
9 to  0 + &9  with  the  incident  beam,  would  be  required.  So 
as  not  to  exceed  a given  angular  resolution,  the  last  counter 
must  subtend  an  angle  no  larger  than  b 9,  However,  the  angle 
0 may  be  as  large  as  2it,  as  is  the  case  with  an  annular  counter. 
Hence,  we  chose  one  dimension  of  the  counter  by  requiring  a 
certain  and  then  made  the  other  dimension  as  large  as  pos- 
sible, consistent  with  an  approximation  to  an  annulus. 

All  the  crystals  were  mounted  in  aluminum  frames  for 
support  with  aluminum  foil,  a few  mils  thick,  being  used 
as  a reflector.  Two  RCA  1P21  photomultiplier  tubes  were 
mounted  at  each  end  of  the  holder  looking  in  at  the  crystals. 
The  outputs  of  the  1P21  phototubes  were  connected  in  parallel. 
The  voltages  supplied  to  the  lP21*s  were  adjusted  with  a 
potentiometer  so  that  pulses  of  equal  height  could  be  ob- 
tained from  a particle  at  either  end  of  the  crystal.  To 
insure  detection  of  all  minimum  ionization  particles,  all 
counters  were  operated  in  the  plateau  region  (II4.OO-I7OO 
volts  on  the  photomultiplier). 

D.  Scattering  Stand 

A picture  of  the  stand  is  shown  in  Figure  2.  The  target 
was  mounted  between  two  screws  (fitting  into  recesses  in 
the  side  of  the  target)  in  the  center  of  the  stand  at  beam 
height  above  the  floor.  On  each  side  of  the  screws  holding  the 
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target  was  an  aluminum  holder  rotatable  through  $60  degrees. 

Attached  to  one  was  a long  pointer  that  swung  about  a large 
protractor  one  foot  in  diameter.  Long  rods  fit  into  the 
aluminum  holders.  X3  and  Xl+  were  rigidly  mounted  and  centered 
on  these  long  rods.  The  rod  could  be  pushed  through  the 
holders  and  then  fixed  in  place  with  set  screws.  Thus  X3  and 
XJl  could  very  conveniently  be  moved  to  any  desired  distance 
from  the  target.  The  angle  of  the  crystal  with  the  beam 
direction  could  be  conveniently  changed  by  rotating  the  holders, 
rods,  and  crystals  all  together,  and  fixing  the  holders  in 
place  with  lock  nuts. 

E.  Electronic  Arrangement 

Figure  3 shows  a block  diagram  of  the  electronic  arrange- 
ment used.  Each  pulse  from  the  scintillation  counters  was 
amplified  and  limited  in  a distributed  amplifier  of  rise 
time  .0026  micro-seconds,  and  a gain  of  10.  This  was  then 
fed  to  a fast  double  coincidence  circuit.  Both  6BN6  and 
bridge  circuits  (Figures  l|  and  5)  were  used  interchangeably. 

With  the  6BN6  circuit  resolutions  better  than  four  milli- 
microseconds were  obtained  and  used  in  some  runs.  However, 
for  most  of  the  work  resolutions  about  1.5  x 10“  seconds 
were  found  sufficient  for  the  ’'fast”  coincidence  circuits. 

Fast  double  coincidences  between  Crystals  1 and  2,  1 and 
3,  and  2 and  l;.  were  made  to  give  optimum  overall  time  resolu- 
tion of  the  selected  total  events.  The  outputs  of  the  coin- 
cidence circuits  were  fed  into  cathode  followers  that  drove  125  ohm 
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RG63U  cable  that  carried  the  double  coincidence  pulses  up 
from  the  cyclotron  building  (3OO  feet)  to  the  laboratory 
building.  Here  they  were  further  amplified  and  fed  to  fast 
discriminator  and  pulse  shaping  circuits  (Figure  6).  The  fast 
double  coincidence  circuits  in  the  cyclotron  building  are  sub- 
ject to  feed-through.  That  is,  a single  very  large  pulse,  per- 
haps due  to  a nuclear  explosion  in  the  crystal,  may  give  an 
output  from  the  double  coincidence  circuit.  It  is  possible, 
however,  to  differentiate  thi3  from  a true  double  pulse  since 
in  our  coincidence  circuits  the  pulse  height  from  double  was 
three  to  ten  times  that  of  any  single  pulse  input.  The  dis- 
criminator circuit  (Figure  6)  was  easily  capable  of  this  dis- 
crimination. The  EFP60  uni vibrator  in  the  output  of  the  dis- 
criminator supplied  uniform  output  pulses  8 volts  in  height 
and  O.lp.  sec.  wide.  The  fast  double  coincident  counts  between 
XI  and  X2  were  used  as  a monitor  for  the  experiment.  This  fast 
coincidence  output,  after  passing  through  the  discriminator 
and  pulse  shaper,  was  sent  to  a fast  0,lp,s  decade  scalar,  and 
also  to  a slower  coincidence  circuit  to  be  combined  with  the 
other  fast  double  coincidences  in  a Rossi  type  (O.lp  sec.) 

" slow"  triple  coincidence  circuit.  The  output  of  this  circuit 
consisted  of  over- all  quadruple  coincidences  in  our  four 
crystals.  Switches  in  the  Rossi  circuit  allowed  us  to  look 
at  the  individual  rates  from  any  one,  or  the  coincidence  of 
any  two  of  the  individual  double  coincidence  circuits  so  that 
we  could  check  for  internal  consistency  of  our  results  and 
"trotible  shoot"  more  easily. 
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III.  EXPERIMENTAL  PROCEDURE 
A.  Analysis  of  the  Beam 

Range  curves  were  taken  to  estimate  the  percentage 
of  a mesons  in  the  beam  through  the  1-2  telescope.  These 
range  curves  were  taken  as  integral  curves,  that  is,  the 
number  of  four-fold  coincidences  for  a fixed  number  of  1-2 
coincidences  was  recorded  as  a function  of  the  total  ab- 
sorber thickness  in  front  of  crystal  number  4*  Typical 
range  curves  are  shown  in  Figure  7 for  mesons,  and  in 
Figure  8 for  n+  mesons.  All  the  points  were  taken  for  the 
same  number  of  monitor  counts.  The  statistical  error  in 
the  point  for  0”  absorber  is  of  the  order  of  1 percent. 

For  purposes  of  analysis,  the  results  were  plotted  as 
a differential  curve  showing  the  relative  number  of  mesons 
stopping  in  each  interval  of  absorber.  The  differential 
curve  corresponding  to  Figure  8 is  shown  in  Figure  9»  The 
roughly  flat  section  of  the  differential  curve  from  0" 
copper  absorber  to  1”  copper  absorber  corresponds  mainly 
to  the  x mesons  which  were  lost  as  a result  of  nuclear  ab- 
sorption of  the  mesons  In  the  copper.  The  rise  in  Figure  9 
after  1”  thickness  of  copper  corresponds  to  the  loss  of 
those  mesons  that  were  slowed  down  due  to  ionization  loss, 
and  came  to  the  end  of  their  range  in  the  copper.  Thus,  most 
of  the  mesons  have  a range  of  about  1.25”  equivalent  of  copper, 
after  leaving  counter  2,  with  a KE  = 80.4  Mev^  and  a momentum 
of  170  Mev/c.  This  range  curve  was  taken  for  the  run  with  a 
I/J4”  aluminum  target  and  positive  mesons.  The  value  of  the 
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energy  during  the  negative  meson  runs  was  similar,  as  the  meson 
target  position  was  unchanged.  However,  for  positive  mesons 
on  a 1/2"  aluminum  target  the  value  of  the  kinetic  energy  was  75 
Mev  after  leaving  counter  2.  For  the  l/4"  and  l/2"  thick  targets 
the  average  distance  before  interaction,  one  half  the  sample 
thickness,  would  result  in  a loss  of  1.2  Mev  and  2.L.  Mev 
respectively.  In  Figure  9 not  many  of  the  mesons  had  a range 
greater  than  1 3/8"  copper,  hence  the  large  drop  in  the  graph 
at  this  point. 

As  mentioned  above,  both  the  hole  in  the  shielding  and 
the  deflecting  magnet  were  essentially  momentum  selectors. 

They  provided  a roughly  monochromatic  it  meson  beam.  However, 
electrons  and  p.  mesons  formed  near  the  Be  target  inside  the 
cyclotron  with  the  same  momentum  as  the  it  mesons,  also  emerged 
from  the  channel.  Since  the  p meson  is  lighter  than  the  it, 
a of  the  same  moment-um  as  the  it  has  a greater  range,  equal 
to  1.73  inches  of  Cu.  Thus  the  rise  above  1 3/4"  in  Figure  9 
corresponds  to  p mesons,  of  the  same  momentum  as  the  it's, 
coming  to  the  end  of  their  range.  These  form  about  4 percent 
of  the  main  beam. 

At  2"  of  absorber,  about  2 percent  of  the  original  positive 
beam  was  still  present,  comprised  most  probably  of  positrons 
and  a few  p»s  produced  from  decay  in  flight  of  it's.  The 
positrons  stopped  very  gradually,  with  increasing  thickness 
of  absorber,  because  of  their  multiplication  via  showers. 

The  electrons  formed  a much  larger  proportion  of  the  negative 
beam.  At  2"  absorber  in  the  negative  range  curves,  about 
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10  percent  of  the  original  beam  remained.  We  assume  that 
this  was  also  the  proportion  of  electrons  in  the  incident 
beam.  The  possibility  existed  that  the  original  number  of 
electrons  in  the  incident  beam  might  be  much  higher,  but 
that  after  penetrating  1.5”  Cu  (about  2.6  radiation  lengths), 
only  a fraction  of  the  original  number  would  be  present.  To 
check  this,  we  examined  Wilson’s  shower  curves  for  200  Mev/c 
electron  in  lead  to  obtain  an  idea  of  the  order  of  magnitude 
of  such  a correction.  Although  we  used  copper  rather  than 
lead,  a rough  comparison  of  the  two  cases  is  possible,  by 
measuring  thicknesses  in  radiation  lengths.  It  was  found 
that,  on  the  average,  each  electron  in  the  primary  beam  pro- 
duced several  electrons  of  sufficient  energy  to  trigger  our 
last  counter  after  traversing  the  absorber.  Further,  in 
one  hundred  cases  examined,  less  than  10  percent  resulted 
in  no  electrons  (i.e.,  all  the  original  energy  being  in 
if-rays)  passing  through  the  absorber.  Hence  this  error  can 
result  in  only  about  a one  percent  error  in  our  final  beam 
composition,  stated  as  approximately  10  percent  electron 
contamination. 

The  range  curves  were  also  influenced  by  the  geometric 
arrangement  of  the  crystals  and  copper  absorber.  Crystal  3 
and  Crystal  1}.  were  placed  apart  with  the  absorber  starting 
at  a distance  3/V  In  front  of  the  last  crystal.  As  copper 
absorber  was  added,  some  mesons  near  the  upper  and  lower 
edges  of  Crystal  3 which  would  have  passed  through  Crystal  1+ 
in  the  absence  of  copper  absorber,  could  now  be  multiply 


22 


scattered  in  the  copper  and  thus  miss  the  last  crystal*  As 
this  effect  was  more  important  near  the  end  of  the  range  of 
the  mesons,  some  of  the  slope  of  the  range  curve  just  before 
the  break  is  attributed  to  this  factor* 

To  better  estimate  the  number  of  mesons  stopping  per 
unit  thickness  of  copper,  we  can  multiply  the  integral  range 
curve  by  a term  e**  to  correct  for  the  nuclear  interaction 
of  the  mesons  in  the  copper.  This,  of  course,  overestimates 
the  beam  left  after  the  71*3  are  stopped  since  the  p-mesons 
have  no  nuclear  attenuation.  Figures  10  and  11  show  these 
plots.  Also  shown  (Figure  12)  is  a range  curve  of  the  ir 
mesons  scattered  by  1/2"  aluminum  at  an  angle  of  35  degrees. 

One  should  note  that  in  the  35  degree  curve  there  is  an  ad- 
ditional energy  loss  of  about  7*5  Mev  in  the  sample.  This 
explains  the  apparent  different  value  of  the  range  of  the  itfs. 
Also,  since  p”  - mesons  and  electrons  were  absent  from  the 
scattered  curve,  the  break  in  the  integral  curve  was  sharper. 
Figure  13  is  a similar  curve  of  ir  on  aluminum  at  160  degrees. 

B * Stability  and  Efficiency  Checks 

During  the  data  runs  frequent  checks  were  made  to  make 
sure  that  the  electronics  were  stable  and  that  the  electronic 
efficiency  had  not  varied.  Thus,  after  a run  at  an  angle  $ , 
the  turret  was  returned  to  the  zero  angle  position  in  the  main 
beam  to  see  if  the  number  of  four-fold  coincidences  for  a given 
number  of  1-2  monitor  counts  had  changed.  This  was  generally 
reproducible  within  statistics  and  provided  a check  both  on 
the  electronics  and  on  the  turret  mechanism. 
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The  electronic  efficiency  wee  measured  by  placing 
Crystals  3 and  I).  in  front  of  Crystal  2,  which  was  smaller  in 
size  than  3 and  4*  Crystal  2 was  also  sometimes  placed  be- 
tween 3 and  Jj..  The  ratio  of  four-fold  coincidences  to  two- 
fold coincidences  gave  the  electronic  efficiency  for  low 
beam  intensities.  At  higher  beam  intensities  the  apparent 
efficiency  decreased,  since  only  the  1-2  circuit  incorporated 
a fast  scalar  (.1  jis  resolution),  while  a slow  scalar  ( Lp  ^ls 
resolution)  was  used  on  the  four-fold  coincidences.  However, 
after  correcting  for  this  resolution  difference  (and  taking 
into  account  the  cyclotron  duty  cycle),  the  results  were 
consistent.  At  other  angles  to  the  beam  direction  the  four- 
fold rate  was  low  enough  not  to  require  correction.  Efficiencies 
might  vary  on  different  runs  due  to  cable  length  variations 
which  were  only  "roughly”  adjusted  to  peak  the  counting  rates, 
but  efficiencies  were  measured  during  each  run,  and  were  con- 
stant for  a given  set  of  cable  lengths  used  during  a run, 
which  might  last  several  weeks.  In  all  cases,  the  efficiencies 
were  of  the  order  of  90  percent  or  greater. 

C • Spurious  Counts.  Accidentals  and  Scattering  In-Out  Rate. 

At  all  angles,  runs  were  made  for  both  target  in  and 
target  cut.  The  results  were  then  subtracted  to  give  the 
number  of  scatterings  due  to  the  target.  With  the  target 
removed,  a few  four-fold  coincidences  appeared  from  several 
sources.  First  there  was  scattering  (nuclear  and  Coulomb) 
from  X2  which  was  placed  close  to  the  target.  In  addition, 
there  was  general  background  which  very  occasionally  gave  a 
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simultaneous  pulse  in  each  of  the  four  crystals.  However, 
a much  more  likely  event  was  the  occurrence  of  background 
oulses  simultaneously  in  X3  and  Xij.  in  coincidence  with  a 
beam  particle  in  crystals  1-2,  since  the  1-2  Instantaneous 
rate  was  as  high  as  1,800,000  mesons  per  minute  (90,000  per 
minute  time  averaged)  during  cyclotron  bursts.  Since  the 
background  intensity  was  higher  near  the  median  plane,  the 
number  of  accidental  coincidence  counts  was  reduced  as  the 
turret  holding  X3  and  XI4.  was  rotated  away  from  the  median 
plane.  This  is  also  why  the  particular  choice  for  fast 
double  coincidence  of  X1-X2,  X1-X3,  and  X2-X1+  was  made.  X3 
and  Xlj.  were  close  to  each  other  (1|.M),  so  if  a stray  particle 
actually  penetrated  X3  and  Xij.,  it  would  have  a high  probability 
of  an  accidental  coincidence  with  the  high  rate  in  XI  and  X2. 

By  pairing  XI  and  X3  and  also  pairing  X2  and  Xlp,  this  possibility 
wa3  reduced.  For  now  the  stray  particle  through  XI  and  X2 
had  to  coincide  within  the  fast  coincidence  resolving  time  with 
the  particle  in  X3  and  XI4.  rather  than  within  the  slow  resolv- 
ing time  of  the  Rossi  circuit. 

Near  the  median  plane,  the  counters  began  to  be  in  line 
with  the  direct  beam.  Thus,  near  0 degrees  angle,  the  four- 
fold rate  actually  decreased  when  the  target  was  inserted. 

This  resulted  from  the  fact  that  the  beam  was  very  well  col- 
limated. When  the  target  was  inserted,  multiple  Coulomb 
scattering  in  the  target  spread  the  beam  out  so  that  it  was 
not  as  intense  in  the  forward  direction.  This  sort  of  effect- 
prevented  us  from  measuring  the  scattering  at  angles  below 
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about  10  degrees.  Here  the  ratio  of  counts  for  "target  in" 
over  counts  for  "target  out"  was  about  2:1.  This  was  due 
to  beam  spread  and  also  due  to  n-mesons  coming  from  the  decay 
of  it  mesons  in  flight  (the  p.*s  from  our  energy  it  meson  had 
a cut-off  angle  of  about  12  l/2  degrees  in  the  laboratory 
system).  At  15  degrees  the  ratio  of  counts  for  "target  in" 
to  "target  out"  was  about  8:1.  For  this,  and  other  reasons, 
we  do  not  trust  the  results  below  15  degrees.  At  the  other 
angles  the  in-out  ratio  varies  from  as  low  as  3:1  UP  to  12:1. 
At  160  degrees  we  had  a ratio  of  about  3.5:1*  During  the 
"target  out"  runs  sufficient  absorber  was  added  between 
counters  3 and  4 to  correct  for  the  increase  in  energy  of  the 
measured  particle  with  the  target  absent.  Table  1 shows  the 
counting  rates  with  3/4-"  copper  and  l/2"  copper  absorber  in 
place . 

D.  Criterion  for  Elastic  Scattering 

In  order  to  measure  the  elastic  scattering,  a copper 
absorber  was  placed  between  X3  and  X4.  Low-energy  protons 
and  mesons  scattered  with  loss  of  energy  would  not  have  had 
sufficient  range  to  penetrate  the  absorber.  The  incident 
beam  was  only  approximately  monochromatic.  The  beam  energy 
spread  on  leaving  counter  2 was  80+5  Mev  (except  for  the 
iC- meson  case  using  l/2"  Al  where  it  was  75+5  Mev).  Hence, 
to  detect  all  the  elastically  scattered  beam  particles  in 
Crystal  4,  the  absorber  could  not  be  too  thick.  The  amount 
of  absorber  chosen  was  such  as  to  stop  mesons  of  less  than 
58  Mev  energy  after  leaving  the  sample  (the  order  of  7*5  Mev 
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fable  1.  H/pical  quadruple  count  imp  rates  as  f unction  of 
absorber*  thickness  for  •sr*  beam,  1/2"  Al.  target. 
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is  lost  traversing  1/2"  of  aluminum) . This  allowed,  of  course, 
inelastically  scattered  mesons  with  a^E-415  + 5 Mev  to  be 
detected  and  counted  as  elastically  scattered  particles.  If 
this  were  an  appreciably  large  part  of  the  cross  section  it 
would  obscure  some  of  the  results,  particularly  in  the  back- 
ward direction  where  the  cross  section  for  elastic  scattering 
dropped  to  its  lowest  value,  and  would  have  to  be  considered 
in  the  interpretation  of  the  results.  Photographic  emulsion 
experiments  indicate  that  inelastic  scattering  with  an  energy 
loss  of  20  Mev  or  less  is  quite  small  But  to  check  this 
result  ourselves,  rough  range  curves  at  a few  angles  were 
taken.  The  35  degree  curve  is  shown  in  Figure  12.  At  small 
absorber  thickness,  inelastically  scattered  mesons,  plus 
reaction  products  such  as  protons,  can  penetrate  the  absorber 
to  give  four-fold  coincidences.  This  number  should  decrease 
as  the  absorber  is  increased.  At  the  scattering  angle  of  35 
degrees  the  uncorrected  range  curve  is  fairly  flat  between 
5/8"  and  7/8"  copper  thickness,  indicating  that  there  is  not 
much  inelastic  scattering  of  n-mesons  in  this  region.  The 
curve  falls  rapidly  at  a range  corresponding  to  that  of  the 
primary  beam.  In  fact,  the  fall-off  is  sharper  and  clearer 
than  the  range  curve  in  the  incident  beam.  This  is  to  be 
expected  as  both  electrons  and  ^.-mesons  in  the  main  beam 
would  not  be  expected  to  scatter  into  large  angles.  Thus  the 
appearance  of  the  plateau  assures  that  at  least  within  the 
spread  of  the  beam  energy  we  do  not  have  appreciable  contribu- 
tions from  inelastic  scatterings  of  more  than  a few  Mev. 
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5.  Additional  Cbeck3 

A farther  check  on  the  nature  of  the  scattered  particles 
was  frequently  made  by  inserting  1 l/2,!  of  copper  between  X3 
and  Xlp  with  the  target  in,  and  1 3/V'  of  copper  with  the  target 
out*  As  can  be  seen  from  the  incident  beam  range  curve,  this 
was  enough  copper  to  cut  out  most  of  the  incident  mesons.  It 
was  to  be  expected,  therefore,  that  the  ratio  of  scattered 
mesons  at  a given  angle  with  3/V*  copper  between  X3  and  Xlj.  to 
the  number  ivith  1 l/2n  copper,  should  be  at  least  as  good  as 
this  ratio  in  the  incident  beam.  This  was  in  fact  the  case, 
thus  further  assuring  us  that  there  was  no  large  spurious  back- 
ground causing  the  counts  measured. 

The  angle  of  scattering  was  measured  by  reading  a pro- 
tractor one  foot  in  diameter  mounted  on  our  scattering  stand. 

The  zero  angle  was  usually  checked  at  the  beginning  and  end  of 
the  runs  to  make  sure  it  had  not  shifted  due  to  slippage  of 
the  pointer  arm  or  some  other  cause.  Actually,  the  stand 
proved  quite  stable,  the  zero  being  reproducible  to  about  an 
eighth  of  a degree.  To  measure  the  zero  angle,  the  number  of 
four-fold  coincidences  for  a given  number  of  1-2* s was  measured 
for  tho  different  angles  through  which  the  turret  was  swung, 

A typical  curve  for  ir  mesons  is  shown  in  Figure  llj.  from  which 
it  is  seen  that  the  zero  corresponds  to  an  angle  of  1 3/k  degrees. 
This  curve  also  gives  an  indication  of  the  incident  beam  width 
which  is  seen  to  be  about  3 1/2  degrees  full  width  at  half 
maximum,  after  which  it  falls  sharply. 
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When  the  target  was  in  place  the  incident  beam  was  further 
spread  out  because  of  the  multiple  Coulomb  scattering  in  the 
target.  The  zero  angle  run  was  then  repeated.  A typical  result 
is  shown  in  Figure  14  for  1/4"  aluminum  and  in  Figure  1$  for 
1/2'*  aluminum.  It  can  be  seen  that  the  presence  of  the  target 
spread  the  beam  out  to  almost  6 degrees  full  width  at  half 
maximum  for  a l/2"  aluminum  target.  This  beam  spread,  curve  can 
be  used  in  making  resolution  corrections  as  it  is  an  empirical 
measure  of  the  beam  spread  including  Coulomb  effects,  TC-meson 
decay,  and  incident  beam  width.  While  most  of  the  runs  were 
taken  with  l/2"  target  thickness,  some  data  points  (with 
poorer  statistics)  -were  taken  with  a l/4"  target  in  an  attempt 
to  reduce  the  beam  width.  The  corresponding  beam  spread 
curves  with  a 1/4"  target  are  shown  in  Figure  14. 

IV.  RESULTS 

A.  Calculation  of  Cross  Sections 

The  cross  section  is  defined  as  the  number  of  interactions 
per  nucleus  in  the  target  divided  by  the  number  of  incident  par- 
ticles per  unit  area.  The  differential  cross  section  in  the 
laboratory  is  defined  as 

do  _ E 

HIT  “oTTO/ 

where  E equals  the  number  of  events  scattered  at  an  angle  9 

n equals  the  number  of  scattering  centers  per  unit  area 
and  is  given  by  Avogadro’s  number  times  the  density  in 
gms./cc  times  the  thickness  of  the  target  in  cm., 
divided  by  the  gm.  atomic  weight  of  the  target  - 
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N equals  the  number  of  incident  particles  - 
^w  equals  the  solid  angle  subtended  by  the  detecting 
counter  telescope. 

It  was  necessary  to  correct  the  raw  values  of  the  above 
quantities  for  Coulomb  scattering;  nuclear  absorption  in  the 
copper  absorber  and  in  the  crystals;  resolution  of  the  scalars; 
for  the  angle  of  the  target  with  the  main  beam;  and  for  electronic 
efficiency  as  discussed  below.  Figure  16  and  Table  2 compare 
the  corrected  results  for  negative  and  positive  mesons  on  1/2” 
aluminum.  The  error?  shown  are  the  standard  deviations  of  the 
number  of  counts  as  it  is  believed  that  the  statistical  errors 
overshadow  the  systematic  ones.  The  angular  resolution  varied 
from  point  to  point  as  we  moved  the  detecting  telescope  in 
from  over  40n  at  the  small  angles  down  to  a minimum  of  about 
10”  at  some  of  the  backward  angles  for  which  there  were  very 
low  counting  rates.  Thus  our  last  crystal  subtended  scattering 
angles  of  + 1.5  degrees  up  to  + 6 degrees.  This  angular  width 
was  superimposed  on  the  intrinsic  beam  spread  with  the  target 
in  place  (as  shown  in  Figure  14)*  The  Intrinsic  beam  width 
was  due  to  the  spread  of  the  primary  beam  and,  in  addition, 
to  the  multiple  Coulomb  scattering  in  the  target.  Thus  the 
Coulomb  scattering  increased  the  full  width  at  half  maximum 
from  three  to  six  degrees  with  a half  Inch  aluminum  target. 

In  order  to  check  the  Coulomb  Interference  region  of  the 
scattering  of  the  positive  and  negative  mesons,  it  was  thought 
desirable  to  reduce  the  incoming  beam  width  and  keep  the  detect- 
ing telescope  further  out  from  the  target.  Hence  runs  were 


EXPERIMENTAL  CROSS  SECTION 
y"  Aluminum  Target 


Table  2.  Experimental  results,  do  for  l/2n  Al.  Target 
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made  with  a l/4”  aluminum  target  with  the  detecting  telescope 
held  at  a distance  of  40M  from  the  target,  for  the  15  to  35 
degree  region.  These  results  are  listed  in  Table  3 and  are 
plotted  in  Figure  17,  With  this  target  the  full  width  at  half 
maximum  of  the  incident  beam  was  4 1/2  degrees  and  the  detecting 
telescope  subtended  an  angle  of  1,9  degrees. 

For  both  targets  used,  points  were  taken  down  to  ten 
degrees.  Points  below  15  degrees  are  not  plotted,  however, 
as  they  are  not  considered  reliable  because  of  the  large  cor- 
rections (which  radically  vary  with  small  changes  in  cyclotron 
operation)  that  would  be  necessary.  These  corrections  are  due 
to  x-jj,  decays  (which  have  for  our  energy  x*s  cut-off  angle  in 
the  laboratory  of  12  l/2  degrees)  resulting  from  the  large 
number  of  n-mesons  in  the  main  beam.  Coulomb  scattering  cor- 
rections and  the  presence  of  main  beam  particles  at  low  angles 
preclude  the  use  of  points  below  15  degrees.  Similarly,  the 
15  degree  point  may  be  considered  less  reliable  than  those  at 
higher  angles, 

B.  Corrections 

The  principle  corrections  were  made  to  allow  for  the  fol- 
lowing: 

(1)  The  counting  of  particles  was  not  100  percent, 
efficient,  since  many  peaked  very  fast  coincidence  circuits  were 
used.  The  efficiency  is  about  90  percent  that  a meson  that 
traverses  counters  1,  2,  3»  and  4 will  give  a four-fold  coin- 
cidence if  it  gives  a 1,2  coincidence.  Thus  tho  observed  1,2, 3, 4 


kz 


counting  rates  must  be  increased  by  about  10  percent,  as 
discussed  on  pages  22  and  27. 

(2)  The  scattered  count  used  in  calculating  cross 
sections  was  the  difference  between  the  scatterer  in  and 
scatterer  out  count  for  3/4"  and  1"  copper  absorber  respect- 
ively, between  counters  3 and-  4«  Range  curves  taken  of  the 
scattered  particles,  at  not  too  small  angles,  showed  that 
there  were  essentially  no  scattered  particles  able  to  traverse 
1 1/2”  copper  (see  Table  l).  Thus,  of  the  total  number  of 
particles  traversing  counters  1 and  2,  only  that  fraction  of 
the  particles  which  can  traverse  1”  copper  absorber,  but  not 
traverse  1 1/2”  copper  absorber,  should  be  included.  S?nce 
the  ^.-mesons  present  were  mainly  of  the  same  momentum  as  the 
it-mesons,  their  mean  range  was  1 3/4”  copper,  (The  electrons 
had  even  greater  ranges  and  are  discussed  separately  later.) 
Thus  an  empirical  check  is  obtained  of  the  fraction  of  the 
1-2  counts  which  should  be  considered  valid  "effective” 
incident  it  meson  flux  by  taking  a range  curve  in  the  forward 
direction  of  the  1,  2,  3,  4 rate  for  a fixed  number  of  1-2' s 
against  copper  absorber  thickness.  If  the  relative  1,  2,  3,  4 
rate  for  no  absorber  is  taken  as  unity,  then  the  difference  in 
1,  2,  3,  4 rate  for  1”  copper  and  for  1 1/2"  copper  gives  the 
fraction  f of  the  1,  2 count  which  should  be  considered  ef- 
fective incident  '^correct  range"x-meson  flux.  Note  that  this 
automatically  takes  account  of  shorter  range  particles  in  the 
incident  beam  due  to  absorption,  inelastic  scattering,  and 
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large  angle  coherent  scattering  in  the  copper  absorber  used 
between  counters  3 and  4 when  measuring  the  scattered  inten- 
sity. The  value  of  f differed  slightly  for  different  runs, 
but  was  equal  to  about  .69  for  the  runs  and  for  the  re- 
runs . 

(3)  The  next  correction  is  due  to  the  net  effective 
angular  resolution  function  in  the  experiment,  including  con- 
tributions due  to  multiple  scattering  in  the  target.  The  net 
total  resolution  function,  including  multiple  scattering  angu- 
lar smearing,  is  measured  directly  by  measuring  the  angular 
spread  of  the  beam  (quadruple  counting  rate  vs  3 "4  telescope 
angle)  in  the  forward  direction  with  the  sample  in  the  beam. 
(Figure  13).  We  need  merely  note  that  any  larger  angle  single 
scattering  event  is  superimposed  on  this  distribxition  since 
the  relative  probability  of  all  other  smearing  effects  is  not 
changed  by  the  presence  or  absence  of  the  large  angle  single 
scattering.  Thus  at  angles  appreciably  larger  than  the 
measured  main  angular  spread  of  the  incident  beam  (after 
traversing  the  sample)  the  experimentally  observed  scattering 
distribution  may  be  regarded  as  the  true  single  scattering 
distribution  viewed  with  a net  experimental  angular  resolution 
function  measured  in  the  manner  described  above.  This  net 
angular  resolution,  which  is  directly  measured  experimentally, 
includes  such  effects  as  the  angular  spread  of  the  beam  on 
leaving  counter  2;  the  widths  of  the  counters;  multiple  scat- 
tering in  the  sample  and  counter  2;  decay  of  the  7t-mesons 
after  scattering  in  the  sample  (the  decay  ti-meson  must  be  able 
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to  traverse  the  copper  absorber  between  3 and  4 to  register, 
so  a modified  x-p,  angular  decay  function  applies). 

The  net  angular  resolution  function  has  the  usual  effects 
in  distorting  the  "true''  ^°/dQ  curve.  If  a Taylor  series  ex- 
pansion of  ^°/dfi  « g(©)  is  made  about  ©=©^ 

g(©,  + x)  « g(©,  ) + g«(®,  )X  + gW(°l)  x2,tt 

l l i 2, 

then  a symmetric  resolution  function  will  give  a measured 
g(©^),  denoted  G(Q^),  which  does  not  depend  on  terms  contain- 
ing  odd  powers  of  x,  but  only  on  g(©^),  and  the  x , x^,  etc, 
terms.  If  g(0)  has  a sharp  dip  at  0^,  averaging  over  a region 
about  ©=©^  will  broaden  the  dip  and  make  it  shallower.  If  g(0) 
is  concave  upwards  at  0=©^,  (*(0^  will  be  larger  than  g(©1). 
Since  the  true  curves  are  probably  concave  upwards  (on  a linear 
plot)  over  most  of  their  range,  Gf©^  will  usually  be  a little 
larger  than  g(©^),  The  curvature  increases  rapidly  at  smaller 
angles  and  the  correction  also  Increases  rapidly.  For  the 
x*  scattering  curve  in  the  region  of  the  interference  dip  be- 
tween nuclear  and  Coulomb  effects  (near  20  degrees)  the  cor- 
rection may  be  particularly  large.  We  present  the  results 
without  correcting  for  this  effect,  however,  since  the  magni- 
tude of  the  correction  is  quite  sensitive  to  the  assumed  shape 
of  G(0),  By  assuming  a hypothetical  G(0)  given  by  the  upper 
curve  shown  in  Figure  18,  and  making  the  necessary  resolution 
correction,  we  obtained  for  g(0)  the  lower  curve  in  Figure  18. 
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V.  DISCUSSION 
A.  "Optical  Model" 

An  attempt  was  made  to  match  the  experimental  results 
using  an  "optical  model".  The  term  "optical  model"  is  often 
employed  to  designate  the  particular  approximation  method 
(of  solving  the  more  general  "optical  model"  problem)  which 
was  employed  by  Fernbach,  Serber,  and  Taylor^  (F.  S.  T.), 

By  "optical  model"  we  here  mean  that,  in  the  most  general 
sense,  the  wave  equation  to  be  solved  for  an  incident  plane 
wave  plus  outgoing  scattered  wave  can  be  represented  by  a 
wave  equation: 

=o 

where  k(r)  is  complex  and  varies  with  position.  We  emphasize 
this  distinction  so  the  approximations  involved  in  assuming 
that  +--he  actual  scattering  can  be  given,  in  principle,  by  a 
solution  of  such  a wave  equation  with  a properly  chosen  k(r) 
will  be  separated  conceptually  from  the  consequences  of 
further  approximations  made  in  the  solution  of  the  problem. 

Such  further  assumptions  could  include  (l)  the  assumption  of 
a uniform  nuclear  model,  i,e,  a constant  nuclear  density 
within  a sphere  of  radius  R,  and  zero  density  outside; (2) 
the  assumption  of  a particular  value  for  R such  as 
R = 1,4  A1^  x lO"1^  cm,;  (3)  the  approximation  of  using  con- 
stant complex  k "inside  the  nucleus",  i.e,  for  r<R,  and  using 
the  k appropriate  for  a Coulomb  field  for  r>R;  (4)  then,  finally, 
one  may  use  various  approximate  methods  of  solving  the  resulting 
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wave  equation,  such  as  the  F.  S.  T.  method,  or  by  use  of  a 
Born  approximation,  etc. 

In  the  case  of  aluminum,  the  above  choice  of  R gives 
kQr  = 3.6,  where  kQ  = k(«>  ).  Thus  the  nuclear  size  is  the 
order  of  a wave  length  rather  than  large  as  compared  to  a 
wave  length.  Assumptions  (2)  and  (3),  although  useful  for 
making  the  problem  definite  and  not  overly  difficult  of  sol- 
ution, are  probably  poor  approximations  to  the  true  situation. 
Since  we  have  no  trustworthy  method  for  estimating  the  addition- 
al distorting  effects  which  are  introduced  by  the  approximation 
methods  (lj.)  above,  we  decided  to  carry  through  exact  phase 
shift  calculations  of  the  definite  model  resulting  from  ap- 
proximations (1),  (2),  (3)  above.  These  calculations  were 
carried  through  for  16  different  choices  of  complex  k » k^  + 
ik£  for  R = 1.1+16  x lO-"^  cm.,  for  both  x+  and  it  mesons, 

k-^  and  were  related  to  an  assumed  complex  potential  + iVg 
using  the  Klein  Gordon  equation 

(E-V)2  = p2c^  +•  -h2c2k2 

2 

where  p =5  meson  rest  mass  and  E = meson  total  energy  = pc  + 
kinetic  energy  at  large  distances.  The  term  quadratic  in 
V was  ignored  for  r>R  so  that  tabulated  Coulomb  wave  functions 
could  be  used.  Values  of  V^,  V^,  k^,  and  k^  used  for  r<R  are 
listed  in  Table  i+.  The  results  of  these  calculations  are 
shown  in  Part  2,  A brief  comparison  of  the  experimental  and 
calculated  curves  shows  that: 

(a)  The  calculated  curves  all  show  strong  diffraction  minima 
in  the  region  of  60  degrees  to  75.  degrees.  The  experimental 
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Table  4.  List  of  Parameters  for  Cases  Studied 


V1 

v2 

kl 

k2 

Case  j'f 

Mev 

Mev 

x 1012 

x 1012 

x 10“13 

1 

-30 

-15 

1 0 J'.O 

0.920 

5.5 

2 

-20 

-10 

9.77 

0.627 

8.0 

3 

-20 

- 5 

9.77 

0.313 

16.0 

4 

-20 

- l 

9.77 

0.0627 

70.0 

5 

-10 

-10 

9.14 

0.642 

7.8 

6 

-5.8 

0 

8.86 

0 

- 

7 

0 

-20 

8.52 

1.317 

3.8 

8 

0 

-10 

8.49 

0.661 

7.6 

9 

0 

- 5 

8.48 

O.33O 

15.2 

10 

+20 

-10 

7.12 

0.714 

7.0 

11 

-20 

-25 

9.81 

1.58 

3.21 

12 

-20 

-35 

9.85 

2.18 

2.29 

13 

-30 

-22 

10.41 

1.35 

3.70 

i'+ 

-30 

-30 

10.44 

I.83 

2.73 

15 

-LlO 

-15 

11.00 

0.904 

5.53 

16 

-40 

-30 

11.04 

1.80 

2.78 
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curves  appear  to  be  much  too  smooth  in  this  region  to  be  con- 
sistent with  this  feature  of  the  theoretical  curves.  The  ex- 
perimental angular  resolution  in  this  angular  region  was  + 
about  7 degrees,  which  will  reduce  the  effect  of  a sharp  dip 
considerably,  but  dips  as  strong  as  the  theoretical  ones  should 
still  show  up  strongly  in  the  experimental  curves.  Considering 
the  statistical  uncertainty  in  the  measured  cross  sections  and 
the  spacing  of  the  points,  however,  it  cannot  be  said  that  much 
less  pronounced  dips  are  certainly  absent.  The  absence  of 

strong  diffraction  dips  is  in  agreement  with  the  present  inter- 

12 

pretation  of  electron  scattering  results  of  Hof stadter,etal. 
and  is  expected  for  a non-uniform  nuclear  model  where  the 
nuclear  density  is  higher  in  the  center  and  tails  off  gradually 
in  an  extended  surface  region.  One  expects  minima  to  reappear 
when  probe  particles  are  used  which  have  a much  shorter  mean 
free  path  for  incoherent  processes  in  nuclear  matter.  Thus  fast 
protons  of  20-100  Mev  would  be  expected  to  show  minima,  and 
n-nesons  should  show  minima  for  larger  nucleii,  where  the 
outer  regions  can  be  more  effective  in  shielding  the  central 
regions.  Schiff^  has  investigated  various  nuclear  distribu- 
tion shapes  in  Born  approximation  to  see  how  gradual  a drop- 
off is  needed  to  remove  the  minima.  The  arguments  are  as  fol- 
lows. If  f&(®)  is  the  basic  scattering  amplitude  for  a 
point  nucleus,  then  fa(9)f^(0)  is  the  Born  approximation 
scattering  amplitude  from  an  extended  nucleus,  where  for  a 
spherically  symmetric  distribution  p(r)  of  nuclear  density 

V®)  * j’”v 
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where 


q = 2k  sin  -g 


is  proportional  to  the  momentum  change  in  scattering,  and 
we  have  tbe  integral  of  ^(r)  normalized  to  unity  over  the 
nuclear  volume.  If  x = qr,  then 


v ■$-  ft* 


If  we  plot  sin  x against  x,  and  superimpose  a plot  of 
x p (q)»  we  can  see  how  inteSranrt  will  vary  with  x.  For 
a uniform  nucleus  f (^)  is  constant  for  x<qR,  and  zero  for 
x>qR.  This  is  shown  in  Figure  19  for  several  choices  of  in- 
creasing q.  Note  that  for  small  q (case  1)  X|*  has  its  entire 
weighting  for  sin  x positive  so  the  integrand  is  always  posi- 
tive, For  larger  q (case  3)  part  of  x^  comes  in  the  negative 
half  cycle  of  sin  x,  and  the  integrand  is  large  and  negative 
near  the  cut-off  value  x^  = qR,  The  net  amplitude  from  the 
integral  will  thus  change  sign  for  some  value  of  x^  between 
x and  2n  and  again  for  x^  between  2n  and  371,  The  diffraction 
minima  occur  when  the  integral  passes  through  zero.  For  the 
uniform  nucleus  the  condition  Is  tan  x^  = x.^  which  has  roots 
for  x-^  = [j.,)|-9,  7«73,  etc.  In  order  to  avoid  having  the  integral 
change  sign,  the  nuclear  density  ^ must  decrease  gradually  for 
large  r so  the  contribution  from  the  regions  x^  = it  to  2it,  etc., 
of  negative  integrand  can  never  exceed  the  contribution  from 

the  regions  = 0 to  it,  etc.  of  positive  integrand.  Results 

12 

of  Hofstadter  etal'  from  electron  scattering  suggest  that  this 
may  be  the  actual  case,  although  the  conclttsions  are  not  cer- 
tain. ^ The  early  interpretation  of  Hof stadter' s results  was 
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based  on  a Dorn  approximation  analysis  arid  favored  a proton 

density  which  decreased  exponentially  with  distance  from  the 

center  of  the  nucleus.  Yennie  etal  have  performed  exact 

phase  shift  calculations  for  a uniform  nuclear  model  using  a 

range  of  Z values.  They  find  that  the  diffraction  minima  are 

similar  to  the  Born  results  for  low  Z but  are  almost  completely 

removed  for  high  Z.  At  the  time  this  paper  was  written,  their 

1 8 

favored  interpretation  for  Z*80  favors  a uniform  nuclear 
model  with  a tail.  The  electron  scattering  for  lower  Z should 
be  more  sensitive  to  these  matters. 

We  note  that  there  is  no  attenuation  shielding  of  the 
central  region  by  the  outer  regions  in  the  case  of  electrons. 

As  the  mean  free  path,  X*  , for  absorption  (incoherent  proces- 
ses) in  nuclear  matter  is  gradually  decreased,  the  inner  region 
of  the  nucleus  will  be  more  and  more  shielded  by  the  outer 
regions,  and  will  have  correspondingly  less  effect  on  the  scat- 
tering. In  a rough  way,  we  may  treat  this  by  saying  that 
p(r)  in  the  expression  for  f^C©)  should  be  reduced  in  amplitude 
progressively  for  decreasing  r by  the  shielding  effect  of  the 
nuclear  matter  at  larger  radii.  In  the  limit  the  effective 
p(r)  differs  from  zero  only  near  the  outer  edge  of  the  nuclear 
distribution,  say  near  rm.  In  this  limiting  case  f-^  changes 
sign  when  xm  = krm  passes  through  x.  This  node  at  3^  = n when 
compared  with  the  unshielded  uniform  nucleus  value  x^  = shows 

that  one  effect  of  the  shielding  is  to  move  the  minima  to  smal- 
ler angles  and  make  the  nucleus  appear  effectively  larger  ac- 
cording to  an  interpretation  which  neglects  shielding  effects. 
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In  Part  2 the  selected  values  of  X^were  all  too  large 

to  have  much  effect  on  the  angle  of  the  minima,  but  it  was 

shtwn  that  fair  agreement  with  the  position  of  the  minima 

using  the  phase  shift  analysis  was  obtained  if  the  Born  analy- 

9 

sis  used  q = 2k-^  sin  ^ , with  the  value  k^  of  the  real  part 
of  the  inside  k,  in  place  of  the  outside  kQ. 

B . Comparison  of  Experimental  and  Calculated  Curves 

(1)  Comparing  the  calculated  n+  and  ir  curves  in  Part  2, 
it  is  seen  that  the  n+  and  ir  curves  for  a corresponding  set 
of  inside  nuclear  potentials  are  nearly  identical  for  9 = 30° 
to  5>0°.  (By  "corresponding”  potentials  is  meant  = -30  Mev, 
for  ir  and  = - 20  Mev,  for  n+,  for  example,  where  a 10  Mev 
difference  in  inside  Coulomb  potential  is  taken  into  account.) 
For  9 <30°  the  n+  curves  show  a destructive  interference  between 
Coulomb  and  nuclear  effects,  while  the  ir  curves  have  a less 
obvious  constructive  interference  effect.  The  principle  dif- 
ference between  the  calculated  curves  and  the  experimental 
curves  in  this  region  is  that  the  experimental  curves  place 
the  it+  and  ir  interference  separation  at  appreciably  larger 
angles  than  do  the  theoretical  curves.  Thus  the  experimental 
separation  extends  from  15°  to  k$° > while  the  theoretical 
curves  place  the  separation  at  about  one-half  or  two-thirds 
of  this  angular  region.  Thus  no  good  quantitative  match  can 
be  made  with  the  it+  dip  using  any  of  the  theoretical  curves. 

For  this  reason,  the  best  choice  of  V for  the  fit  is  made 
using  the  ir  curves.  A detailed  comparison  of  the  experimental 
and  theoretical  ir  curves  shows  that  a fairly  good  quantitative 
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fit  can  be  obtained  over  the  angular  region  20°  - 45°  using 
= -30  Mev  or  = -4-0  Mev  with  the  choice  of  Vg  relatively 
unimportant . 

If  one  observes  that  the  experimental  7C+  and  7 r curves 
match  well  for  9>4-5°  the  difference  at  smaller  angles  may  be 
interpreted  as  a shift  in  the  angular  position  of  the  inter- 
ference minimum  relative  to  those  of  the  theoretical  curves. 

The  71+  curve  could  then  be  said  to  roughly  match  the  choice 
= -20  Mev  or  -3O  Mev.  Increasing  tends  to  damp  the  amount 
of  the  dip  in  the  calculated  ti+  curves  near  15  degrees  and  one 
might  hope  to  select  a favored  value  of  on  the  basis  of 

this  damping.  This  selection  is  rendered  quite  difficult,  how- 
ever, in  view  of  the  great  difference  in  position  of  the  inter- 
ference in  the  experimental  and  theoretical  curves.  The  experi- 
mental curves  of  Figure  17  for  the  l/h"  aluminum  sample  have 
much  better  angular  resolution  than  those  of  Figure  16  in  the 
region  below  30  degrees,  and  a much  better  match  with  the 
theoretical  curves  is  possible.  In  making  the  comparison, 
however,  we  note  that  all  of  the  experimental  points  should 
be  lowered  somewhat.  In  particular,  those  at  1 $ degrees  are 
probably  much  too  high  due  to  contributions  from  still  smaller 
angles  because  of  the  resolution  width.  The  theoretical  curves 
are  shown  corresponding  to  (V.^,  V^)  = (-30,  -22)  and  (-4.0,  -If?) 
f or  tt,  and  (-20,  -10),  (-20,  -2£),  and  (-30,  -If?)  for  7i+.  Com- 
parison shows  that  = -30  to  -34.  (Mev)  is  required  for  the 
7i'  curve.  The  7t+  curve  cannot  be  matched  perfectly  due  to  the 
fact  that  the  theoretical  curves  maintain  low  cross  sections 
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to  much  lower  angles  than  the  experimental  curve.  Reference 
to  Figure  17  and  Figures  23  and  2 ij.,  and  noting  something  like 
the  "area"  of  the  dip,  as  well  as  the  match  for  ©-18°.  sug- 
gests again  the  = -20  Mev  is  best.  The  theoretical  curve 
for  V2  = -25  Mev  seems  overdamped  and  that  for  = -10  Mev 
seems  underdamped,  so  a value  for  between  -10  Mev  and  -25 
Mev  seems  best. 

(2)  It  is  of  interest  that  the  experimental  scattering 
cross  section,  in  addition  to  showing  no  nuclear  size  dif- 
fraction minima,  maintains  a fairly  high  value  near  90  degrees, 
and  actually  rises  some  in  the  backward  direction.  The  optical 
model  phase  shift  and  Born  approximation  calculations  both 

j _ 

give  a steady  decrease  in  ^ between  interference  minima.  We 
note  that  a Born  approximation  calculation  of  the  optical  model 
is  equivalent  to  assuming  isotropic  scattering  from  each 
element  of  nuclear  matter,  so  the  net  f ( Q ) due  to  the  nuclear 
effects  alone  has  just  the  shape  of  the  nuclear  distribution 
form  factor  fb(9).  Another  approach  would  be  to  use  the  fa(9) 
indicated  by  the  meson-nucleon  (x+  + p and  or  + p)  coherent 
scattering.  The  experimental  curves  of  for  + p are  peaked 
in  the  backward  direction  and  the  value  at  180  degrees  is  a 
rapidly  increasing  function  of  the  relative  energy.  We  note 
that,  for  the  backward  scattering  of  a meson  by  a nucleon  in 
the  nucleus  to  be  elastic  and  coherent  with  respect  to  the 
nucleus  as  a whole,  the  nucleon  must,  after  absorbing  the  full 
momentum  transfer,  still  be  in  the  same  quantum  state.  The 
fact  that  this  becomes  increasingly  difficult  as  © is  increased 
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is  just  reflected  in  the  decrease  in  the  calculated  coherent 
form  factor  and  hence  is  automatically  taken  into  account  if 
for  the  elementary  scattering  is  nearly  independent  of  the 
energy  of  the  relative  motion.  The  collisions  which  contribute 
nainly  to  the  coherent  back  scattering  correspond  to  selecting 
Fourier  momentum  components  from  the  bound  nucleon  wave  function 
which  are  directed  towards  the  meson  initially,  and  which  are 
directed  away  from  the  meson  after  the  collision.  Since  the 
nucleon  wave  function  is  not  changed  for  coherent  scattering, 
over- all  energy  and  momentum  balance  is  made,  noting  that  the 
nucleus  aosorbs  the  recoil  as  a unit.  The  feature  to  notice 
is  that  the  main  contribution  to  the  coherent  back  scattering 
comes  from  meson-nucleon  collisions  where  the  relative  motion 
energy  is  much  larger  than  if  the  nucleon  were  stationary, 

j 

Thus  backward  meson-nucleon  — curves  for  110  Mev  or  13£  Mev 
may  be  more  suitable  to  use  than  those  for  80  Mev,  However, 

-3  _ 

there  is  a factor  of  3 to  5 increase  in  ^ at  180  degrees  in 
going  from  78  Mev  to  the  higher  energies.  This  is  capable  of 
giving  a qualitative  explanation  of  the  rise  in  ^ at  larger 
angles,  but  it  should  be  emphasized  that  this  explanation  lies 
outside  of  the  concept  of  an  optical  model,  i.e.,  of  the  use 
of  a wave  equation  with  a complex  k(r;.  It  does  point  a 
direction  in  which  one  might  modify  the  optical  model  results 
by  taking  account  of  the  angular  dependence  f (9)  of  the  elem- 

E 

entary  nucleon  scattering  process.  Multiple  scattering  effects 
will  probably  tend  to  make  some  f ( © ) intermediate  between  that 
for  the  elementary  processes,  and  an  isotropic  distribution 
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corresponding  to  f(0)  be  the  effective  f (Q), 

The  value  of  the  potential  V in  nuclear  matter  can  be 
related  to  the  scattering  amplitudes  in  the  forward  direc- 
tion of  the  elastic  scattering  of  pions  on  protons.  Suppose 
particles  with  propagation  constant  k are  incident  on  a slab 
of  material  with  N scattering  centers  per  unit  volume. 

The  effect  of  the  scattering  centers  may  be  considered 
equivalent  to  the  presence  of  a potential  V inside  the  medium 
and  a propagation  constant  k^  = kQ  + k^"  for  the  particles  in- 
side the  medium.  Consideration  of  the  resulting  wave  function 
for  elastically  scattered  particles  in  the  forward  direction 
after  traversal  of  the  medium  leads  to  the  results 

k2  = k2  + k N f (0 ) 

where  f(0)  is  the  real  part  of  the  coherent  scattering  ampli- 
tude in  the  forward  direction  of  the  incident  particle  from 
one  scattering  center.  Taking  the  Klein-Gordon  equation  for 
the  particle  inside  and  outside  the  medium, 

(E-V)2  = EQ2  + E2c2k12 

K2  = E2  + h2c2k2 

and  combining  with  the  above  equation  for  k^,  we  get 

v = It*  N l»2cg  f (°) 

v v - 

In  the  case  of  meson  scattering  in  nuclei  we  have  not 
one  type  of  scattering  center  but  two,  viz.,  neutrons  and 
protons.  The  scattering  of  it"  and  n+ mesons  off  protons  has 
been  measured  by  Fermi  etal^  Bodansky  etal^  and  using  a 
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charge  independent  isotopic  spin  formalism  (which  is  equiv- 
alent to  assuming  that,  except  for  Coulomb  effects,  the  scat- 
tering amplitude  of  it-  mesons  on  neutrons  equals  the  scattering 
amplitude  of  mesons  on  protons)  Fermi  calculates  the  phase 
shifts  from  which  the  scattering  amplitude  for  elastic  scat- 
tering in  the  forward  direction  may  be  derived. 

Suppose  n mesons  are  incident  on  the  nucleus.  Then  we 
can  calculate  f(0)  for  the  or  on  the  protons.  Since  there 
are  neutrons  present  also,  we  need  f(0)  for  rr  on  neutrons. 

But  we  assume  this  equals  the  amplitudes  from  m+  scattering 
from  protons.  Hence  we  can  consider  in  aluminum  the  elemen- 
tary scattering  centers,  N per  unit  volume,  to  be  neutron 
proton  pairs  with  a scattering  amplitude  equal  to  the  sum  of 
the  scattering  amplitude  for  r on  p and  it*  on  p.  Transform- 
ing the  scattering  amplitude  to  their  values  in  the  laboratory 
system  and  evaluating  V at  65,  120,  and  135  Mev  using  the  real 
part  of  the  scattering  amplitude  to  calculate  the  real  part  of 
the  well  yields  the  curve  shown  in  Figure  20,  It  is  not  com- 
pletely clear  what  energy  should  be  used  in  evaluating  V,  In 
the  meson-nuclear  scattering  experiments  the  laboratory  energy 
is  used  and  there  is  no  net  momentum  transfer  for  forward  scat- 
tering, For  a bound  nucleon,  there  are  several  pos  sible  com- 
plicating effects  of  the  binding.  The  fact  that  the  nucleons 
are  moving  would  be  unimportant  if  the  basic  scattering  ampli- 
tudes were  slowly  varying  with  respect  to  the  meson  energy  in 
a system  where  the  nucleon  is  at  rest.  This  is  not  true,  how- 
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ever,  and  the  nucleon  momentum  distribution  should  be  con- 
sidered. For  a nucleon  moving  toward  or  away  from  the  meson 
with  20  Mev  kinetic  energy  and  t^I.02)  an  incident 

meson  of  momentum  P and  total  energy  E will  have  momentum  and 
energy  P'  and  E’  in  the  Lorentz  system  where  the  nucleon  is 
at  rest,  where  cp>  = *n(cp  +/&«,E)  and  E*  = *n(E  + {*n cp).  In 
this  case,  rnBnE  4-0  Mev*  Thus  some  sort  of  weighted  average 
of  V over  a region  2+0  Mev  to  120  Mev  is  required.  Since  the 
curve  of  V vs  E is  roughly  linear,  this  correction  should  be 
relatively  small  and  an  effective  E of  80  or  90  Mev  is  prob- 
ably reasonable.  The  kinetic  energy  inside  the  nucleus  prob- 
ably should  be  increased  by  the  amount  of  the  effective  change 
in  potential  energy  on  entering  the  nucleus,  ivhich  would  favor 
a value  of  20  to  30  Mev  higher.  If  this  is  done,  values  of 
45  Mev  are  obtained. 

A second  unknown  effect  of  the  binding  is  the  possibility 
of  changing  the  effective  reduced  mass  of  the  meson-nucleon 
system  and  thus  appreciably  changing  the  elementary  scattering 
strength.  This  effect  is  almost  always  ignored  and  the  con- 
ceptual features  of  the  possible  effect  are  often  confused 
with  an  effective  reduced  mass  effect  on  the  amount  of  phase 
space  available  for  scattering.  To  clarify  the  concept,  at 
least,  let  us  consider  the  scattering  of  slow  neutrons  by 
protons  bound  in  a chemical  lattice.  For  the  singlet  inter- 
action, the  relative  wave  function  for  a square  well  poten- 
tial behaves  as  sin  K r inside  the  range  of  the  potential  and 
does  not  reach  90  degrees  of  phase,  so  the  scattering  length 
is  negative.  In  considering  the  scattering  of  bound  protons. 
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it  Is  assumed  that  there  Is  no  change,  due  to  the  binding, 

-12 

in  the  behavior  of  the  relative  wave  function  for  R<^10  ,e.g. 

If,  however,  the  proton  were  bound  by  a spring  of  such  stiffness 
that  a motion  of  10“^  Cm.  gave  a LV  of  ^10^e#v,  there  would 
be  a significant  increase  in  the  reduced  mass  of  the  relative 
motion,  and  the  value  of  K inside  the  range  of  the  force  would 
be  correspondingly  increased.  This  could  lead  to  a bound 
singlet  state  of  the  deuteron,  and  to  a positive  scattering 

amplitude.  For  atoms  in  a chemical  lattice,  the  bV  is ^ 1 e.v 

-9 

for  "'lO  ' cm.  motion,  so  such  effects  are  negligible.  For 
nucleons  in  nuclear  matter,  however,  it  is  not  clear  that 
such  effects  could  not  occur.  If  they  did  occur,  the  net 
effect  might  still  be  small  unless  a near  resonance  condition 
existed,  as  for  the  singlet  system.  It  is  not  yet  clear  whether 
meson-nucleon  interactions  should  be  considered  as  having 
resonances,  or  whether  the  concept  of  an  effective  interaction 
potential  defining  the  scattering  is  appropriate.  It  is_  ap- 
propriate to  note  that  such  effects  are  not  easily  excluded 
on  general  grounds  and  should  be  remembered  as  an  interesting 
possibility. 

Referring  to  Figure  20,  it  is  seen  that  the  value  of  V 
for  E = 80  Mev  is  3^  Mev,  which  is  only  approximately  in  agree- 
ment with  the  value  selected  previously  in  matching  the 
theoretical  and  experimental  scattering  curves. 

Inside  the  nucleus  we  have  a complex  propagation  constant, 

+ iKg,  -^e  wave  Unction  inside  the  nucleus  is  pro- 
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portional  to  e lr  so  that  the  complex  propagation  constant 
introduces  an  attenuation  of  the  incident  wave  function  pro- 
portional to  e"^2r#  The  beam  intensity  is  proportional  to 
the  wave  function  squared  so  that  the  beam  will  be  attenuated 
by  a factor  e”2^2r  corresponding  to  a mean  free  path  for  ab- 
sorption b\  g A.  table  of  ^>s  for  each  of  the  cases 

calculated  in  Part  2 is  shown  in  Table  ij.  together  with  the 
absorption  cross  sections  calculated  using  the  phase  shift 
analysis.  The  calculated  absorption  cross  section  as  a func- 
tion of  the  mean  free  path  for  absorption  is  plotted  in  Figure 
21.  The  calculated  points  are  shown  as  dots  and  some  of  the 
curves  drawn  for  different  values  of  have  been  extrapolated 
in  the  region  of  long  mean  free  path.  If  we  choose,  as  above, 
the  value  V-^  = 30  Mev  for  ai”mesons  and  V'1  = 20  Mev  for 

17 

mesons,  by  comparing  with  measured  absorption  cross  sections  1 
(which  may  have  to  be  corrected  for  large  angle  elastic  scat- 
tering) we  can  choose  a mean  free  path  A0*  By  considering 
the  published  absorption  cross  sections,  and  correcting  for 
large  angle  elastic  scattering  where  necessary,  we  choose  a 
best  present  value  of  530  millibarns  for  the  absorption  cross 
section  of  ir  on  aluminum,  and  we  choose  a 10  percent  smaller 
value,  or  lj.80  millibarns  for  the  it+  value.  The  upper  solid 
horizontal  line  in  Figure  21  is  drawn  at  the  value,  the 
lower  at  the  x+  value.  The  mean  free  path  for  absorption  for 
both  negative  and  positive  mesons  is  then  seen  to  be  about 
5 x 10“^cm.  Note  that  a 10  percent  change  in  the  absorption 
cross  section  would  change  the  mean  free  path  by  a bit  over 


64 


1 x 10_^cm.  It  should  be  kept  in  mind  that  the  mean  free 
path  calculated  by  this  technique  is  of  course  only  as  good 
as  approximations  1-3  above.  Thus,  for  example,  if  the 
nuclear  radius  were  smaller,  then  the  measured  absorption 
cross  sections  would  be  closer  to  geometric  cross  section 
and  yield  a smaller  mean  free  path.  Above,  by  matching  the 
theoretical  curves  to  the  experimental  curves  in  the  region 
of  small  angles  we  saw  that  a value  for  equal  to  or  less 
than  25  Mev  (corresponding  to  X = 3*2  x 10-^cm)  would  fit 

a 

the  data  but  that  a value  of  — 10  Mev  (corresponding  to 
X = 8 x 10”^3  cm)  was  definitely  too  low.  This  second 

9. 

method  of  determination  of  the  mean  free  path  thus  indicates 

a value  of  X between  about  3 x 10~^  cm  and  6 x 10**^^. 

8. 

Xa  can  also  be  estimated  from  the  elementary  pion  nucleon 
scattering  data  in  two  different  ways.  First,  we  can  cal- 
culate the  imaginary  part  of  the  potential  inside  the  nucleus, 
exactly  as  we  did  the  real  part,  V^,  on  page  57  hy  using 
the  imaginary  part  of  the  forward  scattering  amplitude.  This 
results  in  values  for  Vg  of  -10  Mev  and  -37  Mev  at  incident 
tc  energies  of  65  Mev  and  120  Mev.  Then  at  80  Mev  incident 
energy  we  would  have  V = (37  + il8)  Mev  while  at  100  Mev 
incident  energy  V = (I4.3  + i27.j?)  Mev.  These  correspond  to 
about  a mean  free  path  of  ij.,6  x 10“^cm.  and  3 x 10~^cm. 
respectively  (see  Table  ij.) . These  values  would  be  lowered  if 
effects  due  to  absorption,  i.e.,  disappearance  of  the  meson, 
are  included. 

We  can  also  consider  the  absorption  coefficient  as  in 
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P.  S.  T. : 


That  is,  the  absorption  coefficient  in  nuclear  matter  equals 
the  particle  density  times  the  cross  section  for  scattering 
of  the  ti  meson  by  a particle  in  the  nucleus.  This  calcula- 
tion gives  a mean  free  path  at  65,  80,  and  120  Mev  of 
5.5  x 10“^cm.,  ij.#2  x 10“^-^cm.,  and  1,9  x 10”^cm.  respectively. 
Hence  at  80  and  100  Mev  incident  energy  we  would  obtain  mean 
free  paths  of  ij.,2  x 10“^^cm.  and  about  3 x 10“^cm.  However, 
the  cross  section  for  inelastic  processes  for  bound  nucleons 
in  a nucleus  will  be  smaller  than  the  total  cross  section  for 
free  nucleons  due  to  the  Pauli  exclusion  principle  require- 
ment that  the  final  nucleon  state  be  previously  empty.  Thus 
these  last  values  should  be  somewhat  increased.  A best 
speculation  as  to  the  true  mean  free  path  might  be  taken  as 
L\.  x 10  ^cm.  for  all  the  above  calculations  with  the  particular 
choice  of  nuclear  model  and  nuclear  size. 
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PART  II 

Calculations  of  the  Expected  Elastic  Scattering  of  ft  Mesons 

by  Nuclei 

Recently  the  present  author  and  others  have  undertaken 
an  experimental  study  of  the  elastic  scattering  of  ti*  and  li- 
me sons  by  nuclei.  At  present,  the  main  investigations  have 
used  aluminum  and  lithium  as  the  sample  materials,  and  80 
and  130  Mev  mesons.  The  experimental  studies  were  under- 
taken with  considerably  greater  angular  resolution  and 
statistical  accuracy  than  earlier  measurements  of  this  type^,^,^‘. 
Since  measurements  of  this  type  can,  in  principle,  yield 
considerable  valuable  information  on  the  interaction  of 
pions  with  nuclear  matter,  we  believed  that  it  was  important 
to  compare  the  experimental  results  with  as  exact  theoretical 
predictions  as  would  be  feasible.  This  has  led  us  to  an 
examination  of  the  methods  currently  favored  for  the  calcu- 
lation of  such  scattering.  We  have  also  carried  through  an 
extensive  program  of  calculations  of  the  scattering  expected 
on  the  basis  of  an  "optical  model'1  using  an  exact  phase  shift 
analysis  for  various  complex  indices  of  refraction  inside  the 
nucleus,  and  using  Coulomb  wave  functions  outside  the  nucleus. 
These  calculations,  involving  many  computer  hours,  provide 
numerous  test  cases  to  compare  with  the  results  of  approximate 
methods  usually  employed  for  this  purpose.  Several  significant 
features  of  disagreement  were  found  between  the  results  of  our 
exact  calculations  and  the  usual  approximate  methods  which  we 
believe  should  be  emphasized  as  giving  important  limitations 
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’ on  the  applicability  of  the  approximate  analysis.  These 

matters  are  discussed  in  some  detail  in  the  following 
sections . 

I 

A second  general  method  of  approach  to  the  elastic  scat- 
tering of  fast  particles  by  nuclei  is  based  on  the  Born  ap- 

i 

proximation,  either  directly  or  through  the  use  of  various 
modifications  which  are  intended  to  include  effects  which 
are  present  in  the  scattering,  but  are  not  readily  incorp- 
orated in  the  usual  first  order  Born  approximation  treat- 
ment. As  is  discussed  in  more  detail  below,  the  formulae 
resulting  from  this  second  approach  include  many  features 
. not  present  in  the  optical  model  approach,  but  which  we  be- 

lieve should  be  present  to  some  extent  in  a more  exact  theory 

of  elastic  scattering.  Similarly,  the  concepts  of  a change 

* 

in  phase  relations  due  to  the  real  part  of  the  index  of  re- 
fraction inside  the  nucleus  differing  from  the  outside  value, 
and  the  concepts  of  extinction  effects  due  to  a non-infinite 
mean  free  path  for  absorption  (non-coherent ) processes  in 
nuclear  matter  are  special  features  of  the  optical  model 
which  are  absent  from  the  simple  Born  approximation  analysis, 
but  would  be  expected  to  be  present  to  some  extent  in  a more 

■ exact  theory.  We  have  tested  various  recipes  for  altering 

■ the  usual  Born  approximation  formulae  to  consider  these  lat- 
ter effects  and,  by  treating  cases  which  were  also  treated 
exactly  by  the  phase  shift  analysis,  are  able  to  make  detailed 
comparison  with  the  exact  results  from  the  optical  model  phase 
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shift  analysis.  These  calculations  and  comparisons  are 

presented  in  detail  in  the  following  sections.  We  shall 

not,  in  this  paper,  be  mainly  concerned  with  the  funda- 
20 

mental  theory  used  to  establish  the  approximate  validity 
of  the  differential  equation  which  forms  the  basis  of  the 
optical  model.  Rather  we  shall  emphasize  the  comparison  of 
the  predictions  of  the  various  methods  of  analysis  discussed 
above  using  what  may  be  considered  a semi-empirical  approach. 

The  Optical  Model 

The  introduction  of  the  use  of  an  optical  model  to  the 
investigation  of  the  elastic  nuclear  scattering  of  fast 
particles  was  mainly  due  to  Pernbach,  Serber,  and  Taylor. ^ 

They  showed  that  fast  neutron  elastic,  interaction,  and  total 
cross  sections  when  analyzed  on  this  basis  seemed  to  provide 
an  excellent  basis  for  interpretation  of  the  experimental 
results.  In  carrying  through  their  derivation  of  actual 
final  formulae  for  comparison  with  experiment,  however,  they 
made  two  important  simplifications  which  are  themselves  not 
implied  by  the  concept  "optical  model",  but  which  facilitated 
the  analysis.  These  were:  (1)  A "uniform"  nuclear  model  was 

assumed  as  has  been  customary  in  most  treatments  of  nuclear 
processes.  This  considers  a spherical  nucleus  of  radius  RQ 
having  constant  nuclear  density  inside  and  zero  density  out- 
side. (2)  An  approximate  method  was  used  to  solve  for  dif- 
fraction effects,  similar  in  philosophy  to  the  usual  approximate 
methods  of  dealing  with  interference  and  diffraction  effects 
in  physical  optics,  as  opposed  to  exact  solutions  in  terms  of 
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differential  equations  and  boundary  conditions*  By  adopting 
this  approximate  method,  they  were  able  to  explore  the  general 
behavior  of  the  theoretical  predictions  over  an  extended  range 
of  the  free  parameters  to  an  extent  that  would  require  a pro- 
hibitively large  program  of  calculations  by  the  exact  phase 
shift  method.  Thus  they  were  able  to  present  an  excellent 
preliminary  survey  of  the  subject  which  could  gradually  be 
improved  by  testing  it  against  exact  calculations  for 
selected  sets  of  parameters. 

As  a result  of  the  general  excellence  of  their  paper, 

however,  there  has  developed  some  confusion  as  to  what  one 

means  by  the  term  "optical  model" * The  authors  find  that 

many  of  the  physicists  with  whom  they  have  discussed  the 

matter  have  considered  (1)  and  (2)  above  as  implied  features 

of  the  "optical  model".  However,  recent  experimental  studies 

21 

of  mesonic  x-ray  transitions,  and  the  recent  studies  of 

12 

the  elastic  nuclear  scattering  of  fast  electrons,  at  Stan- 
ford in  particular,  indicate  that  the  nucleus  has  a con- 
siderably higher  density  of  nuclear  matter  near  the  center 
than  previously  believed  to  be  the  case,  with  a gradual 
dropping  off  of  the  density  in  the  outer  regions  (fuzzy 
edge).  Since  the  basic  "optical  model"  differential  equa- 
tions are  capable  of  exact  solution  by  a phase  shift  analysis 
of  the  various  angular  momentum  components,  the  method  need 
not,  in  principle  or  fact,  limit  itself  to  the  approximate 
method  of  Fernbach,  Serber,  and  Taylor. 

7«e  define  the  "optical  model"  as  follows:  The  particle 

being,  scattered  is  represented  at  a large  distance  in  terms 
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of  a plane  wave  and  an  outgoing  radial  wave  from  the  scat- 
terer,  which  we  take  as  centered  at  the  origin  of  coordi- 
nates. (When  a Coulomb  field  is  present,  the  plane  wave  is 
modified  in  the  usual  way.)  At  all  points  the  wave  equation 
for  the  particle  can  be  represented  as 

V 2y  + (kx  + ik2)2Y  = 0 (1) 

where  k^  and  are  real  functions.  In  a region  of  zero 
potential  energy  k^  = kQ  and  k^  = 0,  so 

K1  + ik2  = ^o  (2) 

introduces  a complex  index  of  refraction  n.  If  absorption 
(incoherent)  processes  do  not  occur  at  r,  and  the  momentum 
is  real  (positive  kinetic  energy),  then  k2  = 0,  and  k^  dif- 
fers from  kQ  due  to  a change  in  momentum  (scalar  value) 
from  its  value  at  large  distances.  When  absorption  (inco- 
herent) processes  occur  at  r,  k2  is  positive  and  is  the 
mean  free  path  for  absorption  (of  intensity  rather  than 
amplitude).  When  r represents  a classically  disallowed 
point  (negative  kinetic  energy  or  imaginary  momentum)  the 
situation  becomes  more  complex.  When  the  non-relativistic 
Schroainger  equation  is  used,  k^  and  k2  are  related  to  the 
real  and  imaginary  parts  (V^  + iV2)  of  the  potential  at  r 
and  to  the  energy  E. 

+ ik2)2/2m  = (E-V^-iVg)  (3) 

If  the  Klein-Gordon  equation  is  used,  then 

c2n2(k1  + ik2)2  = (E-V^-iVg)2  - (MC2)2  (4) 
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where  m is  the  rest  mass  and  E the  total  energy.  VJe  have 
used  (i+)  for  relating  k-^  and  kg  to  a (V-^  + iVg)  for  nuclear 
matter  in  treating  the  it  meson  scattering.  VJe  normally  con- 
sider situations  where  the  scattering  system  is  very  massive 
relative  to  the  incident  particle  and  can  thus  consider  lab- 
oratory and  center  of  mass  coordinates  as  equivalent.  Elastic 
coherent  scattering  requires  that  there  be  no  change  in  the 
internal  wave  function  of  the  scatterer. 

In  the  phase  shift  analysis,  ^ is  expanded  in  spherical 
coordinates  yielding  separate  radial  equations  for  each 
angular  momentum  term.  If  ^ corresponds  to  the  part  of  r V 
for  orbital  angular  momentum  1,  and  if  we  neglect  spin  ef- 
fects, the  equation  for  'f  is 


O 


(5) 


with  the  boundary  condition  'f  - 0 at  r ~ 0,  and  Y|  and  *(% 
continue  at  boundaries. 

In  the  case  of  no  Coulomb  potential,  this  solution  at 
large  distances  behaves  as  sin  (kQr  - — ) and  the  scat- 
tering amplitude  is  then 


^ ^ H (C6* 

Vo  |»0 


(6) 


In  the 

pi ex  and  the 

<3\ 


ease  where  there  is  absorption, 
total  non-coherent  cross-section 


will  be  com- 
is 


(7) 


When  a Coulomb  term  is  present  outside  the  nucleus. 
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there  is  Coulomb  scattering  modified  by  the  nuclear  inter” 
action*  The  analysis  is  similar  in  the  relativistic  and 
non- relativistic  cases  if  (a)  terms  quadratic  in  V are 
neglected  in  (I4.)  for  the  region  outside  the  nucleus,  (b) 

p 

if  E/c  is  used  in  place  of  the  rest  mass  m for  the  outside 

Coulomb  wave  functions.  The  theory  in  this  case  can  be  found 
22 


in  Schiff. 


His  equation  20.24  and  20.10  give: 

00  - 


f(Q)  = f (9)  + 


S'" 


1*0 


-\a(  % *»> 


fc(e)  “ 2k  , 20 

o sin  “2~ 

Pin 

where  e p £I+5L  - \ <*) 


(8) 


(9) 


2 

2 1^2  6 /r  the  Couloinb  potential  and  v is  the  velocity 
of  the  particle.  For  a point  Coulomb  field,  the  behave  at 
large  r as  sin  (kQr  - -N; — -«<  log  2 kQ  r + n^  ) showing  the 
logarithmic  variation  of  the  phase  shift  at  large  r due  to 
the  l/r  type  potential.  When  the  extended  actual  nucleus 
is  considered,  an  extra  phase  shift  is  added  at  large  r. 

To  actually  carry  out  the  calculation  of  the  values, 
the  model  of  a uniform  density  nucleus  was  used  with  a con- 
stant complex  potential  + iV ^ inside.  The  logarithmic 
derivatives  of  the  inside  and  outside  solutions  were  matched 
at  r = 1^.  The  inside  radial  solutions  are  of  the  form 
where  ^2*  9 c‘*s  ^ = ^ , 

and  functions  for  higher  1,  and  the  derivatives  of  the  functions 
are  obtained  from  the  recursion  formulae  for  Bessel  functions 
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^4+  and  ^ ^ * Here  ^ s(V»*  1 

and  (k^  + ikgjR  = p.  + iv,  where  we  use  pQ  = kQRo  = 3.600 

for  aluminum  and  79  Mev  x mesons  (kQ  = 8.5* ••x  10^/cra  and 

R = 4.2  x lO^^cm.  The  momentum  was  1?0  Mev/c.).  The 

regular  and  irregular  Coulomb  wave  functions  were  taken  from 

the  tables  of  Bloch,  Hull,  Broyles,  Bouricius,  Freeman,  and 
21 

Breit,  and  the  tables  of  Coulomb  wave  functions  prepared 
by  the  National  Bureau  of  Standards. ^ The  method  of  using 
these  tables  is  not  completely  straightforward,  and  it  has 
been  outlined  in  the  appendix  for  reference. 

Two  methods  were  used  by  Fernbach,  Serber,  and  Taylor. 

In  the  one  usually  employed  for  comparison  with  experiment, 
the  wave  front  is  considered  to  pass  through  the  nucleus  with- 
out disturbing  initial  ray  directions,  and  the  relative 
phase  and  amplitudes  of  points  on  a wave  front  are  considered 
after  traversing  the  nucleus.  All  parts  not  in  the  geometric 
shadow  of  the  nucleus  have  the  same  amplitude  and  phase.  A 
point  on  a part  of  the  surface  for  which  the  ray  traversed 
nuclear  matter  has  an  additional  phase  J (k1-kQ)ds,  and  is 
attenuated  by  a factor  exp  - fkgds  in  amplitude.  (Note; 

Our  kg  is  half. of  their  K,  which  refers  to  intensity  attenu- 
ation.) The  general  diffraction  formula  is  now  used  and  the 
problem  treated  as  the  linear  combination  of  a case  where 
the  amplitude  and  phase  on  the  wave  front  are  the  same  in 
the  geometric  shadow  as  elsewhere,  (giving  no  scattering) 
plus  the  case  where  the  amplitude  is  zero  outside  the  shadow, 
and  equal  to  the  (phase)  vector  change  in  amplitude  due  to  the 
nucleus  inside  the  shadow.  The  scattered  radiation  pattern 
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is  then  that  due  to  a flat  circular  area  and  depends  only 
on  kRQ  sin  9.  In  particular,  it  is  symmetric  about  90  degrees. 

p 

The  absorption  assymptotically  approaches  7iRQ  as  K— 9 . 

It  is  evident  that  the  above  angular  pattern  for  the 

scattering  can  only  be  valid  at  small  angles  since  symmetry 

about  90  degrees  requires  that  all  odd  angular  momentum  terms 

be  £ero  in  (6).  Also  the  Born  approximation  calculation, 

which  should  become  exact  in  the  limit  of  weak  interactions, 

depends  on  the  vector  change  in  momentum  in  scattering, 

q = 2kQ  sin  9/2.  Thus  an  improvement  might  be  expected  In 

their  scattering  formula,  by  replacing  sin  9 with  (2  sin  9/2). 

They  mention  another  formula,  (9),  which  is  based  on  a WKB 

phase  shift  analysis  and  involves  a series  in  P (cos  9). 

This  would  be  expected  to  extrapolate  better  to  large  angles. 

The  errors  in  this  method  were  also  discussed  by  Pasternack 
29 

and  Snyder. 

Recently  their  approximate  expression  for  o^  in  terms 
of  XQ  has  been  used^*^  in  the  analysis  of  meson  scattering 
to  obtain  an  estimate  of  the  mean  free  path  for  absorption 
of  it  mesons  in  nuclear  matter  from  the  experimental  o*  . 
Comparison  of  the  approximate  formula  for  o^with  the  results 
of  our  exact  calculations  shows  good  agreement  for  small 
ok/tiRo  , but  serious  disagreement  for  larger  K,  since  the 

p 

exact  calculation  allows  to  exceed  7tRQ  for  relatively 
small  values  of  KRQ.  The  results  are  compared  in  Tables  5 
and  6.  The  discrepancy  may  be  considered  as  due  to  the 
small  number  of  51  values  required.  InP.  S.  T.  a summation 


i 

O 

r- 1 

0 

T3 

m 

cd  © 

AJ 

XA 

H 

xO 

CA 

O' 

CA 

rH 

Pt 

sO 

sO 

— i. 

XA 

rH 

C 

EH  *H 

• 

• 

0 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

O 

* <M 

— i. 

Al 

pt 

pt 

1 

E — 

XA 

rH 

rH 

O' 

rH 

O' 

XA 

C^- 

CO 

•H 

(h  tH 

7^ 

-a 

o 

-d- 

rH 

pi 

CC 

CA 

Pt 

i-l 

sO 

CA 

CA 

CO 

• 

-P 

© x3 

pt 

CA 

AJ 

CA 

AJ 

rH 

AJ 

XA 

x£ 

XA 

sO 

Pt 

sO 

O 

42  2 

® 

U E 

1 

co 

© 

W 0 

c n 

co 

{ 

o 

o 

u 

© 

/ 

o 

AJ 

sO 

pt 

O' 

O' 

xO 

XA 

rH 

CA 

CO 

AJ 

XA 

O' 

c 

0 

• 

• 

• 

• 

• 

• 

• 

• 

0 

• 

• 

• 

• 

0 

C 

u 

CO 

SO 

Pt 

co 

O 

1 

-d- 

XA 

rH 

O 

CA 

O' 

CO 

sO 

XA 

Pt 

o 

© o 

A- 

•sO 

CA 

co 

rH 

co 

O 

3 

CO 

A- 

pt 

sO 

> 

•H 

j *i 

AJ 

rH 

AJ 

-d“ 

AJ 

rH 

CA 

Pt 

Pt 

CA 

Pt 

1 

O 

1 

cd 

m* 

«s 

i 

© 

• 1C, 

vO 

(-: 

XA 

O' 

O 

rH 

O 

O 

CA 

O 

o 

sO 

rH 

0O 

XA 

XA 

■P 

too  O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

! 

(3 

"EVsl^ 

CO 

rH 

O 

o 

XA 

O 

XA 

CA 

CA 

XA 

rH 

CA 

O' 

XA 

O' 

XA 

AJ 

\ 

•rH 

O' 

sO 

rH 

XA 

AJ 

CA 

O' 

A- 

AJ 

A- 

AJ 

CO 

CA 

CO 

O' 

to 

3-  *« 

Pt 

CA 

AJ 

CA 

-d 

AJ 

i — ! 

AJ 

XA 

sO 

XA 

3 

XA 

xO 

O 

rH 

C 

•H 

V 

© 

1 

P> 

-P 

CO 

t) 

XA 

co 

rH 

rH 

AJ 

XA 

XA 

O 

O 

A- 

• 

CO 

II 

0 

0 

• 

• 

• 

• 

• 

0 

0 

•p 

*X 

CA 

A! 

rH 

O 

AJ 

o 

rH 

CA 

i— ! 

o 

co 

1 

C 

pt 

© 

© 

XA 

A- 

'O 

O' 

CO 

XA 

00 

A- 

Pt 

GO 

AJ 

Pt 

xO 

O 

AJ 

CA 

Pt 

• 

• 

• 

• 

• 

• 

0 

• 

• 

0 

• 

• 

• 

• 

• 

• 

o 

II 

<A 

CO 

O' 

H 

XA 

o 

CA 

xO 

C'- 

CA 

A- 

A- 

o 

Al 

CA 

CA 

o 

o"C 

CA 

1 — 1 

i — ! 

Pt 

XA 

Pt 

o 

Pt 

A- 

O' 

CA 

co 

6 

•d 

( — 1 

sO 

xO 

XA 

XA 

c — 

Pt 

O' 

rH 

A- 

A- 

O 

O 

O' 

E 

CA 

• 

• 

• 

0 

• 

• 

0 

• 

• 

• 

• 

• 

• 

• 

0 

0 

co  || 

AJ 

CO 

CO 

o 

CO 

o 

XA 

Pt 

>o 

XA 

AJ 

xO 

xO 

co 

CA 

AJ 

XA 

•3 

XA 

CO 

_H- 

rH 

xO 

O' 

XA 

AJ 

CA 

xO 

CO 

CO 

o 

O' 

-d- 

O 

/ 

Ih 

cd 

42 

rH 

i — 1 

rH 

rH 

AJ 

rH 

Al 

CA 

O 

-P 

iH 

rH 

vO 

o 

o 

O 

o 

XA 

00 

AJ 

a; 

O' 

xO 

Pt 

O 

XA 

O 

XA 

T3 

rH  CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

0 

0 

© 

•H  II 

— i_ 

o 

AJ 

XA 

o 

O 

Pt 

rH 

O' 

XA 

CA 

A- 

O 

Pt 

OO 

-P 

E _ 

3 

O' 

AJ 

AJ 

XA 

o 

vO 

o 

O 

i— i 

o 

rH 

CC 

rH 

rH 

cd 

H 

H 

rH 

rH 

rH 

AJ 

AJ 

At 

Al 

1 — 1 

Al 

AJ 

r- 1 

£ 

2 

•H 

A 

o 

p 

rH 

O 

XA 

CA 

Pt 

QO 

O' 

XA 

rH 

O' 

Pt 

H 

pt 

Al 

sO 

Pt 

rH 

cd 

0 rH 

• 

0 

• 

• 

• 

0 

0 

0 

• 

• 

• 

• 

• 

• 

• 

0 

O 

II 

CA 

AJ 

rH 

AJ 

O' 

o 

rH 

AJ 

O' 

CO 

O' 

XA 

oo 

A- 

vO 

AJ 

rH 

, 

O' 

ao 

XA 

rH 

00 

AJ 

O' 

XA 

c — 

rH 

AJ 

o 

rH 

CO 

rH 

CA 

£ 

ot 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

• 

XA 

Pt 

CA 

O 

H 

00 

XA 

CA 

-d- 

O 

O 

[ — 

'O 

sO 

O 

A- 

* 

O 

• 

• 

• 

• 

o 

• 

• 

0 

• 

• 

• 

• 

• 

• 

• 

• 

© 

II 

XA 

A- 

PT 

sO 

xO 

a- 

C — 
rH 

pt 

O' 

O' 

A- 

CA 

CA 

rH 

o 

CA 

AJ 

rH 

AJ 

CA 

AJ 

CA 

CA 

CA 

CA 

Pt 

-4. 

aJ 

EH 

AL, 

it 

* 

© 

rH 

© 42 

rH 

Al 

CA 

pt 

XA 

xO 

A- 

00 

O' 

O 

rH 

CM 

CA 

XA 

vO 

+ 

co  E 
cd  3 

rH 

rH 

rH 

H 

rH 

rH 

rH 

U a 


Table  6.  n+ Calculated  Total  Minus  Coherent  Scattering  Interaction  Cross  Sections 


7 6 


CQ 

c 

Ph 

aJ 

o-r' 

•H 

rH 

o 

CO 

O' 

XA 

rH 

O' 

O 

<A 

CO 

rH 

rH 

-d 

a- 

O' 

XA 

i — 1 

0 

o 

• 

• 

• 

• 

• O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

rH 

A- 

-d 

tA 

O 

o 

<A 

CM 

CM 

(A 

C — 

sO 

XA 

CO 

CM 

O' 

CM 

•H 

ii 

CO 

XA 

rH 

XA 

AJ 

CM 

O' 

A- 

A' 

-d 

O' 

nO 

tA 

CM 

-d 

O' 

e 

-d 

tA 

CM 

tA 

-d 

CM 

rH 

CM 

XA 

1A 

1A 

'O 

XA 

sO 

O 

rH 

sO 

c^- 

o 

-d 

vO 

rH 

<A 

O' 

-d 

-d 

vO 

o 

O 

CO 

XA 

tA 

-d 

• 

• 

• 

• 

• O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

II 

~d 

<a 

( — 

rH 

rH 

CO 

-d 

-d 

1A 

CM 

CM 

XA 

XA 

o 

-d 

(A 

oc 

CVJ 

rH 

rH 

rH 

rH 

tA 

-d 

tA 

_d 

(A 

XA 

O' 

(A 

A- 

o 

CVJ 

XA 

XA  <A 

IA 

-d 

o 

O' 

CVJ 

o 

rH 

(A 

O 

O' 

• 

• 

• 

• 

• O • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

<A 

CM 

rH 

— L. 

O' 

rH  XA 

CO 

XA 

r-* 

O' 

rH 

XA 

XA 

nO 

rH 

XA 

II 

-d 

CO 

3 

'O  CO 

-d 

CM 

cA 

-d 

C^- 

A- 

O 

CO 

<A 

O 

Of 

H 

rH 

rH 

rH 

rH 

r-H 

CM 

(A 

CM 

A- 

rH 

a- 

rH 

O' 

O' 

CO 

co 

vO 

o 

O' 

nO 

O 

c— 

sO 

XA 

• 

• 

• 

• 

• o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

II 

O' 

CVJ 

vO 

-d 

CM 

o 

O' 

-d 

o 

o 

-d 

o 

CO 

rH 

CM 

CO 

CO 

XA 

O' 

CVJ 

tA 

sO 

o 

A- 

rH 

r-H 

1 — 1 

1 — 1 

CO 

i — 1 

rH 

rH 

rH 

( — l 

rH 

H 

CM 

CM 

CM 

CM 

rH 

CM 

CVJ 

CO 

CM 

CO 

C^- 

-d 

vO 

A- 

CM 

CO 

XA 

a- 

1A 

XA 

tA 

A- 

rH 

rH 

• 

• 

• 

• 

• o 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

CO 

O' 

CO 

rH 

CO 

CM 

(A 

O 

CM 

sO 

CM 

-d 

CO 

CM 

CO 

rH 

ii 

CO 

A- 

-d" 

rH 

CO 

CM 

O' 

'O 

CO 

rH 

CM 

O 

<A 

CO 

o 

fA 

<5* 

rH 

rH 

rH 

rH 

i — 1 

1 — 1 

rH 

<A 

o 

rH 

CM 

rH 

CM 

CO 

C^- 

XA 

O 

CO 

o 

sO 

rH 

A* 

o 

II 

<5* 

• 

rH 

-d 

• 

CO 

CM 

• 

o- 

rH 

• 

-d 

• 

O 

CJ 

° A 

vO 

tA 

• 

sO 

CM 

• 

vO 

rH 

• 

<A 

tA 

• 

O' 

<A 

• 

o 

-d 

• 

o 

-d 

• 

3 

• 

• 

-d 

• 

tA 

-d 

© © 
M .Q 

d S 

rH  CM  <A  _d 

XA  'O  A-  CO 

o o 

rH 

CVJ 

<A 

-d 

XA 

vO 

o 0 

S3 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

+ 

Vs 


77 


over  ^ values  is  replaced  by  an  integral,  and  this  approximation 
fails  when  only  the  first  few  ^ values  are  involved  in  the 
series. 

In  our  phase  shift  calculations  k R = 3.600,  so  it 
would  be  customary  to  neglect  terms  for  1 ^ Ip.  To  be  certain 
that  no  errors  were  involved,  however,  we  calculated  terms 
through  Q = 5»  and  it  was  evident  from  the  results  that  contri- 
butions from  1-6  would  really  be  negligible.  The  calculated 
angular  distribution  of  the  cross  section  has  been  calculated 
for  the  16  choices  of  complex  potential  each  for  tc+  and  x~  listed 
in  Table  I4..  Tables  of  values  of  the  amplitude  coefficients 
K,  L,  N,  and  M for  selected  angles  (see  equation  A-8)  are 
given  in  Tables  7>  8 and  9 from  which  all  of  these  curves 
can  be  calculated.  Tables  10  and  11  also  yield  the  amplitudes 
and  do/dO  for  the  cases  considered.  It  was  found  to  be 
quite  useful  to  make  vectox-  plots  of  f = jy  vs.  Q of  the 

type  shown  in  Figures  31  to  35  to  interpolate  values  of 
do/dO  between  the  values  calculated.  This  was  particularly 
true  in  the  region  of  the  diffraction  minima,  where  the 
exact  angle  of  the  minimum,  and  the  minimum  cross-section, 
could  be  readily  found.  Also,  it  served  as  a check  on  the 
overall  calculations  since  the  resulting  curves  behaved  in 
a regular  fashion  when  no  errors  were  made.  If  an  occasional 
point  seemed  to  be  out  of  line  with  the  general  curve,  an 
extra  check  was  initiated  and  the  error  found  and  corrected. 

The  results  for  a number  of  cases  are  shown  in  Figures  22  to  30. 


Table  7.  Coulomb  Amplitude  Functions 


2 k f 

c 

= K + i L (ir) 

= -K  + i L M 

6 

K 

L 

10 

14,59880 

-7.29940 

15 

6.80188 

-2.58873 

20 

3.93068 

-1.18705 

25 

2.57326 

-0.62465 

30 

1.81676 

-O.35639 

45 

0.84267 

-0.081702 

60 

0.49575 

-0.015044 

65 

0.42950 

-O.OO53511 

70 

0.34647 

+0.0049134 

75 

0.33449 

+0.0061893 

90 

0.24763 

+O.OI3789 

105 

0.19632 

+0.016565 

120 

0.16440 

+0.017475 

135 

0.14420 

+O.OI768I 

150 

0.13173 

+O.OI7623 

160 

0.12664 

+0.017567 

170 

0.12371 

+0.017518 
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Table  8. 


:ase  # 

Mo 

Mi 

m2 

M 

3 

m4 

m5 

1 

-.72656 

-.76711 

-.75199 

-.24288 

-.02347 

-.001423 

2 

-.44 433 

-.50867 

-.43375 

-.11583 

-.01207 

-.000612 

3 

-.29901 

-.35472 

-.27262 

-.040287 

-.00072 

.000009 

4 

-.15192 

-.21006 

-.11196 

.025729 

.00925 

.001422 

5 

-.36112 

-.43603 

-.32596 

-.10360 

-.01697 

-.002794 

6 

-.007780 

-.010772 

.001115 

.0019382 

.0004 

.000014 

7 

-.61674 

-.76707 

-.53778 

-.19197 

-.0344 

-.001751 

8 

- .44046 

-.52704 

-.33356 

-.11028 

-.02119 

-.004097 

9 

-.29609 

-.35427 

-.20746 

-.064656 

-.01186 

-.001765 

10 

-1.01514 

-1.05142 

-.50212 

-.14008 

-.02221 

.0007 

11 

-.68033 

-.81185 

-.73411 

-.29035 

-.04726 

12 

-.74356 

-.93249 

-.85531 

-.36903 

-.06857 

13 

-.77567 

-.85425 

-.84071 

-.32983 

-.04356 

14 

-.80357 

-.93067 

-.91667 

— .40166 

-.06414 

15 

-.97512 

-.95117 

-I.OIO39 

-.36320 

-.02522 

16 

-.91551 

-1.01162 

-1.05860 

-.50147 

-.07367 

H0 

N1 

N2 

N3 

n4 

n5 

1 

.20641 

.59353 

.55265 

.31449 

.06316 

.006497 

2 

.13843 

.47182 

,448i4 

.19875 

.03645 

.00384 

3 

.22362 

.55120 

.53752 

.20488 

.03298 

.002915 

4 

.30967 

.64347 

.62254 

.20642 

.03189 

.003703 

5 

-.13324 

.15773 

.15702 

.068508 

.01248 

.000163 

6 

-.32894 

-.07486 

.00479 

.OO6033 

.00099 

.000025 

7 

-.27639 

-.03191 

-.079293 

-.026205 

-.00066 

-.000061 

8 

-.40183 

-.17088 

-.10515 

-.030896 

-.00526 

-.001568 

9 

-.4907 

-.26828 

-.13685 

-.038615 

-.0067 

-.001476 

10 

-.6012 

-.52587 

-.45183 

-.16809 

-.02908 

.001542 

11 

-.03504 

.34101 

.23501 

.13343 

.03676 

12 

-.099055 

.30878 

.13684 

.078631 

.03203 

13 

.075173 

.49615 

.40859 

.25292 

.06200 

14 

-.026921 

.42371 

.28518 

.I8386 

.05761 

15 

.23176 

.68575 

.64214 

.45737 

.09549 

16 

-.033214 

.48254 

.34225 

.25158 

.08400 
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Table  9. 


k+ 


Case  # 

Mq  M1 

m2 

m3 

1 

-1.09558  -1.28571 

-1.11599 

-.38645 

-.048827 

-.OOij.288 

2 

-.87256  -1.05011 

-.82857 

-,22?28 

-.029413 

-.002642 

3 

-.85368  -1.03291 

-.79879 

-.17851 

-.021612 

-.001915 

4 

-.83281  -1.03263 

-.76579 

-.13349 

-.01574? 

-.001693 

5 

-.60587  -.75596 

-.51728 

-.14135 

-.020033 

-.001789 

6 

-.24532  -.33858 

-.16643 

-.03383 

-.005032 

-.000496 

7 

-.63958  -.78109 

-.47062 

-.1437 

-.022713 

-.001216 

8 

-.42912  -.52998 

-.31196 

-.08687 

-.013316 

-.001497 

9 

-.28694  -.35218 

-.19940 

-.05212 

-.OO7636 

-.00053 

10 

-.59985  -.53006 

-.16055 

-.02498 

-.00174 

-.000955 

11 

-.87103  -1.11045 

-.84705 

-.30728 

-.048577 

12 

-.85553  -1.15259 

-.84593 

-.32618 

-.057776 

13 

-1.01653  -1.25724 

-I.05309 

-.40905 

-.059173 

14 

-.94899  -1.24289 

-1.00118 

-.41493 

-.068242 

15 

-1.26260  -1.46887 

-1.34489 

-.56226 

-.065686 

16 

-1.00817  -1.33499 

-1.13286 

-.52725 

-.O86683 

No  N1 

N2 

N3 

n4 

N5 

1 

.55139  .38243 

.10681 

.15569 

.033977 

.00323 

2 

.74101  .52150 

.30489 

.1519 

.02524 

.002474 

3 

.90917  .68619 

.49876 

.2009 

.02903 

.002375 

4 

1.0793  .84864 

.69311 

. 24044 

.033403 

.003414 

5 

.64692  .41321 

.20121 

.07009 

.00993 

-.000197 

6 

.80857  .6499 

.36292 

.08996 

.012569 

.001124 

7 

.36009  .034412 

-.12642 

-.04605 

-.005902 

.000228 

8 

.42869  .15842 

.OI678 

-.0003 

-.01652 

-.001584 

9 

.49120  .25317 

.09475 

.02147 

.002308 

-.000393 

10 

-.14183  -.49612 

-.35173 

-.10752 

-.014814 

.002118 

11 

.44776  .20759 

-.077116 

.OOO78 

.007765 

12 

.35934  .092952 

-.23510 

-.09540 

-.005543 

l? 

.43399  .26032 

-.054216 

.05832 

.023384 

14 

.35192  .15947 

-.19104 

-.03480 

.010689 

15 

.36007  .29028 

-.051193 

.I8241 

.050711 

16 

.25634  .12906 

-.26859 

-.05261 

.020486 
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For  comparison  also  carried  through  a number  of 
Born  approximation  type  calculations,  using  various  modifi- 
cations of  the  usual  Born  approximation  procedure  to  try  to 
obtain  a better  agreement  with  various  features  of  the  phase 
shift  calculation  results.  The  various  procedures  are 
listed  Born  approximation  1 to  £,  and  are  discussed  below. 
Born  Approximation  1 (modified) 

We  use  f (0)  for  a point  nucleus 
0. 


fa(e)  = - g^1!2-—  + \ Rc2  tvx  + iv2)^  do) 

> *0 

For  fg(0),  the  nuclear  distribution  form  factor,  use 
q^  = 2k^  sin  0/2  rather  than  qQ  =*  2kQ  sin  0/2.  Here  kQ  and 
k = k^  + ikg  are  the  values  of  k at  r = , and  inside  the 

nucleus  respectively.  The  use  of  k^  rather  than  kQ  in  fg(0) 
brings  the  diffraction  minima  to  about  the  same  angle  as 
the  phase  shift  calculation. 

Also  in  fg,  we  weight  interior  regions  of  the  nucleus 
less  than  the  surface  by  a factor  e 2'  o”  ' to  try  to  take 


account  of  the  attenuation  effects.  Thus 


f = ^ — 
B 


■3  ( 


(11) 


Born  Approximation  2 (standard) 

Here  we  use  equation  10  for  f (0)  and  use  the  regular 

8. 


Born  approximation  also  for  fg(0),  Thus,  this  is  the  true 


first  Born  Approximation. 

r - 3 [ 2 sin  q„r 

j r t — - 


■B 


V 


) dr 


(1?) 
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Born  Approximation  3 (modified) 

Here  use  equation  10  for  fa(d)  and  use  q = 2 (k^  + ikg) 
sin  9/2  for  fg(0),  thus  q is  complex. 


fB  Ro3 


(sin JLS-J  dr 


q r 


(13) 


Born  Approximation  k (modified) 


Here  we  do  not  separate  f(Q)  into  a product  of  f and 

BL 

fB,  since  the  Coulomb  contribution  outside  from  the  central 
protons  is  not  attenuated,  even  if  their  short  range  force 
effect  is  attenuated. 


(a) 

f = 

f + f 
Z1  xo 

due 

& 

(b) 

1! 

•H 

2m  ( 

"P  > 

(c) 

It 

O 

S-l 

2m 

" P* 

L 

( 


o 

sin  q r 

^ or 

q r 
Ho 


(34) 


) (Vx  + iV2)e-k2(Ro“r)dr 


*1* 


1*2 


e2 


dr 


Born  Approximation  5 (modified) 

Here  we  used  the  same  technique  as  in  approximation  ij., 
but  without  the  attenuation.  This  differs  from  the  true 
Born  approximation  (number  2)  in  that  the  potential  is  taken 
as  constant  for  r<RQ,  while  in  number  2,  the  true  inside 
Coulomb  potential  also  appears.  Thus  this  approximation 
matches  the  phase  shift  calculations  in  holding  V fixed  for 
r <Rq. 


The  results  of  the  calculations  are  partly  available 
from  inspection  of  the  Figures  25,  26  and  30,  In  general, 
the  separation  into  f fg  always  gave  a true  zero  at  the  inter- 
ference minima  except  in  approximation  3,  when  a complex  q 
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was  used.  VJhen  th^  outside  qQ  was  used,  the  diffraction 
minima  occurred  at  angles  independent  of  the  choice  of 
and  V^,  while  the  phase  shift  calculations  gave  minima 
at  angles  that  decreased  as  became  more  negative..  The 
angles  of  the  minima  were  well  matched  by  the  Born  approxi- 
mation calculation  using  rather  than  qQ  for  fg. 

Approximation  3 represented  one  attempt  to  obtain  the 
proper  damping  of  the  diffraction  minima.  The  calculations 
showed  considerable  over  damping  in  this  case.  Approximation 
4 and  5 also  should  not  have  true  zeros,  but  they  were  much 
less  damped  than  the  phase  shift  calculations  results. 
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Appendix  Outline  of  the  Procedure  for  Making  the  Phase  Shift 

Calculations 


Equation  5>  can  te  written  for  r > RQ, 


lA-1)  W -o 


Z.  Z?e 

For  Z = kr  and  °<  = — 

R <\r 

The  solutions  of  this  equation  are  given  by  the  hypergeometric 
function.  The  solution  regular  at  r = 0 is  denoted  F9  (•<  ,Z), 
and  the  irregular  solution  is  denoted  by  Gj  (•<  , Z).  These 
solutions  are  tabulated  by  Breit  etal^^  (called  A),  and  in 
"Tables  of  Coulomb  Functions"  published  by  the  National  Bureau 
of  Standards^  (called  B).  The  functions  are  normalized  such 
that  they  have  asymptotic  forms  for  large  z. 


(a)  Fj  (*,  Z)  sin^  Z -M  In  2 Z -*#■  + ">0 
(A-2) 

(b)  G^(<,  Z^-*cos£z  -«tln  2 Z - V-  + Vl 

where  = Arg  V + 1 + i«<  ). 

For  a pure  Coulomb  field  (non-relativistic ) only  F.  can 
appear  to  satisfy  the  boundary  condition  at  r = 0.  For  our 
problem  of  matching  inside  and  outside  solutions  at  r = R, 
we  use  a linear  combination  of  Fj  and  G|  which  leaves  the 
incoming  radial  wave  the  same  as  for  Fq  alone.  This  requires 
that 


(A-3)  \ 

V, 


^F|_  ooaSji  + Gj  sin  £*"]  , or 

F^  + ^ G$  4-  iF|^  e'^  sin  Sjj 


Ill 


where  ^ = additional  nuclear  phase  shift  so 

(a-4)  tf!  -* 

as  Z-*  °°  . Matching  logarithmic  derivatives  of  the  inside 
and  outside  solutions  at  r * RQ  gives 


outside 


(A-S)  + lk2)  **  llllnslde  = ko 

The  inside  solutions  are  of  the  form 

(a)  ^ (Z)  = Jz  J|  + 2/2  where 

(b)  fz 

J-l/2  ^ S C0S  Z 

(A-6) 

(c)  ^ Ji/2  (z^  “ 3in  Z and 

(«)  't/  ‘ 4 

Z = (k.^  + ikg)  r and 
(kx  + ik2)  R 0 = H + i* 

The  left  side  of  equation  A-5  is  evaluated  numerically 
for  \ - 0 to  5 from  the  above  formulae.  To  evaluate  the 
right  side  of  equation  A-5,  we  proceed  as  follows; 

«-tr,  £. 


(A-7)  k 


¥*  - k T *v  + (Gt  + IB*  ) e1**  sin 

o^  ko  LF|  + (gJ  + iljj  ) el^sinSj 

Solving  for  and  substituting  in  equation  8 can  be 
written 

(a_8) 


where  K,  L,  N,  M are  real.  Here  K and  L are  the  real  and 
imaginary  parts  of  ft(9),  for  which  must  be  obtained. 
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Tables  of»<yo  vs,e<  are  given  in  reference  B (note  the 
differences  of  notation  as  shown  in  Table  ll}.,),  Values  of 
K and  L for  selected  angles  are  given  in  Table  7»  Solving 
for  N and  M reduces  to  evaluating  Fj  , F|  , , G^  at 

the  surface • 

To  evaluate  Fj  and  we  follow  the  procedure  of 
Reference  B where 


(A-9) 


(a)  Fj  = C%  ( <*  ) Z*  +1  $ (*  , Z)  where 

/o2  . *2x1/2 

(b)  cx  (*  ) = l] — Cs-1  (<  } 


Tables  of  C ( °*  ) for  positive**  are  given  in  Reference 
B.  For  negative  •< 


(A-10)  C,  (-•*)  = 6**  CA  ( H ) 

The  functions  ^ are  directly  tabulated  in  Reference 
B,  Interpolation  procedures  are  given  to  obtain  accurate 
values  of  ^ for  values  of  W and  Z between  those  tabulated, 

(A-n)  ^(**v.*>  »\.T h <?*> + ' ' ' 

where  tables  of  X*' ^ , Z)  are  given  for  n = 1,2,3,,, 

F^  is  given  by 

* . 0 S ft  Uil  1* 


(a) 

(A-12) 


■S-ov  % V 0 . 


Fov  *s  °>  *** 

>1A 


(b)  IV* 4 F,1  - “ 1 F,t| 
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To  determine  G^  we  use  the  relation 

. „ F,  Gj-i 

(A-13)  *»’-? — IT 

This  requires  knowing  G^  or  Gq.  It  turns  out  to  be 
easisr  to  obtain  G1  from 

(A-14)  = A cos  ^ x 

Reference  A gives  Tables  of  A^  and  ^ (Tables  31  and  36) 
in  terms  of  ci  and  Z. 

To  obtain  G^  , we  use  the  Wronskian  relation  to  solve 
for  G^  in  terms  of  other  quantities 

(A-15)  Fj  Gjj  - F|  = 1 

We  are  now  in  a position  to  evaluate  the  right  side  of 

equation  A-7  in  terms  of  e^*  sin  , which  is  required  in 

2i  i 

equation  8.  Actually  we  found  it  easier  to  solve  for  e^1  * 
which  is  simply  related  to  e^  * sin  # We  also  require 
the  values  which  are  readily  obtained  from«»|0  by  noting 

p (l  + X + W)  = (*♦  ^ (^  ’*0 

so 

V’  Vi  + 

♦ Tables  of  N and  M for  selected  Q for  our  cases  are 

given  in  Tables  8 and  9» 

i 

> 


that 
(A-l6 ) 
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